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Hydrocarbon fementation is a subject of intensive 
research because of its anticipated potential of pioducing 
Single^dell protein. The system consists of hydrocarbons dispersed 
in the aqueous medium in the form of drops, aqueous medium itself, 
microbial cells, and air bubbles. Stirred vessels with specially 
designed agitators for vigorous mixing of the different phases 
are commonly used for this purpose. Hydrocarbon drops continuously 
coalesce and redisperse giving rise to a distribution of drop 
sizes. Simultaneously, the microorganisms consume hydrocarbons 
either directly from the drop-in terf ace or those dissolved in the 
aqueous medium. Therefore a thorough knowledge of drop- inter action 
processes and iiie droi>-size distributions is essential before a 
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rational design and scale-up of laydrocarbon fermentors may be 
attempted. 

The pcesent study has attempted to characterize the nature 
of drop- interactions and the drop-size distributions obtained in 
hydrocarbon fermenters and related systems. A coalescence- 
redispersion (G-R) model has been proposed for drop- interactions 
and its predictions have been verified for a large number of' 
experimental data obtained during the course of this investigatio 
and from the literature. Experimental measurements for drop- 
interaction fhequencies have been carried out using a new method 
based upon the C-R model. The model has finally been used, 
together with the growth equations, for a Monte-Carlo simulation 
of growth of biomass in hydrocarbon fermenters. 

Chapter II describes a method for experimental measurement 
of drop-size distributions in liquid-liquid agitated systems. A 
detergent stabilization method was chosen because of its appli- 
cability to drop-size measurement in hydrocarbon - water systems 
both in the absence and the presence of biological grovi/th. This 
method consists of sampling a small amount of dispersion in a 
mildly agitated concentrated solution of detergent followed by 
micro pictography. A high degree of reproducibility could be 
obtained by this method. All the experimental work was carried 
out in a three-litre fementor fitted with a Waldhof- agitator. 

Drop-size distribution data experimentally observed in 
gas-oil and in n-alkane fermentations were obtained from various 
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sources (Prokop and Ludvik, 1973; Blanch and Einsele, 1973; and 
Ka linger, 1973) and an attempt was made to empiricallj correlate 
them with experimental conditions. It was found that the number 
density distributions for different experimental conditions of 
speeds of rotation, dispersed phase fractions, sizes of vessels, 
and differaat times of fermentation, when normalized with the 
Sauter-mean diameter, reduce to a single normalized distribution. 
Thus the Sauter-mean diameter emerged as a possible basis for 
scale-up of hydrocarbon fermenters. The Sauter-mean diameter was 
correlated with jHrameters describing experimental conditions in 
the circulation- stirring systems. The exponent of dependence 
upon Weber number in particular was found to be -(0,1-0. 3) instead 
of -0.6 obtained for flat-blade turbine systems. Another obser- 
vation which emerged out of the studies was that the. concentration 
as well as the physiological state of the biomass strongly affect 
the Sauter-mean diameter. 

In Chapter III, a coalescence-re dispersion model for drop- 
interactions is discussed. The model views the drop- phenomena as 
a single step coalescence-redispersion process between any two 
droplets which results in redistribution of m_ass between them. 

The population balance equation for batch vessels at equilibrium 
drop-size distribution is amenable to analytical solution and for 
the case of constant coalescence frequency, the model predicts an 
exponential equilibrium distribution of drop sizes. Exiperimen tal 
data for drop-size distributions, collected in labor.atory during 
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the present inTestigation for n- alkane-water system mixed by a 
Waldhof-agitator, and for various other systans available in the 
literature (Brown and Pitt, 1972| Luhning and Sawistowski, I97I1 
Shinnar, 19611 and Prokop and Ludvik, 1973), were satisfactor ily 
predicted by the model. On comparing the model predictions with 
those of the drop-size distributions proposed by vanHeuven and 
Hoevenaar (1968) and Bayens (1967), a sign test rated the C-R 
model better than the latter two. The model predictions for 
Waldhof-agitated systems (or circulation stirring systems) were 
found to improve when the coalescence frequency was assumed to 
be a function of drop- sizes. 

Based upon the C-E model, a method was proposed for the 
measurement of coalescence frequency. It requires measurement of 
transient- state drop-size distributions in a stirred batch vessel, 
or steady state DSDs of in- and out-going streams in a continuous 
flow stirred tank vessel. Prom the data obtained in the batch 
vessel, tbe coalescence frequency can be calculated by the slope 
of the straight line obtained by plotting the second moments of 
drop-size distributions vs time lapsed since the start of the 
stirring. In this connection, the experimental data were obtained 
for n- alkane- water system agitated by a Waldhof-mixer. The results 
show that instead of a power function dependence of coalescence 
frequency upon the speed of rotation as normally observed for 
flat-blade turbine systons, the Waldhof-agitated systems 
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are characterized by a Linear dependence of drop-in teraction 
frequency on the speed of rotation. 

I’inally, the C-R model has been used to. represent the drop- 
phenomena in hydrocarbon fermentors Vv'hile attempting to predict 
their behavior, A Monte-Carlo simulation procedure has been 
utilized for this purpose and the results are presented in 
Chapter IV, The nature and significance of interactions between 
dispersed phase d3x>ps has been investigated. The results of 
computer simulation confirm the empirically observed self preserv- 
ing nature of drnp-size distributions in batch hydrocarbon 
fermentors. A study of the contribution of continuous phase 
growth towards overall biomass growth has also been attempted. 


* * * 



CHAPTER I 


IRTROnJCTIOH 

World-Wide protein shortage has stimulated a vigorous 
search for new sources of protein supply, in the last fifteen 
years. By the end of this century, the demand for protein is 
estimated to he around 65 million tons, of which only 43 million 
tons is expected to he met hy conventional means (Vilenchich 
and Akhtar, 1971). In order to overcome such high deficiencies, 
it is evident that every possible means he tapped to the utmost. 
Microorganisms which contain proteins as much as sixty percent 
of their dry weight, will play an important role in meeting the 
requisite demand in future. A strong case fbr their cultivation 
is presented hy high rates of multiplication of microorganisms, 
relatively small floor-space requirements and a high efficiency 
of conversion of substrate to proteins. In addition to hydro- 
carbons, other carbon substrates for the growth of microbes are 
methanol, ethanol, natural gas, and glucose. Hydrocarbon fermen- 
tation has, however, the additional advantage that it results 
in improved quality of the petroleum product used for the supply 
of C-substrate. In India, there is a potential for producing 
6 lakh (600,000) tons of pure proteins from Assam crudes alone 
(Iyengar, £t a^. , 1968). 

Hydrocarbon fermentation involves four phases which 
contribute to the growth of microorganisms. These are the gas 
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phase (air) , the aqueous phase (taediutn) » the dispersed phase 
(oil),, and the solid phase (cells). Transfer of oxygen and 
O-suhstrate across interfaces, sometimes more than one, before 
being consumed by the cells, makes the analysis of these systems 
extremely complicated. However, we shall restrict our attention 
in the present study to substrate limited growth alone. In other 
words, oxygen concentration in the system is assumed to be above 
the critical level all the time. - Hydrocarbons being sparingly 
soluble in aqueous medium, exist in the form of drops. These 
drops are present in a turbulent environment, and continuously 
coalesce and redisperse. The cells consume the hydrocarbons 
either directly from, the interface or those dissolved in the 
continuous phase:. The substrate limited growth in hydrocarbon 
fermenters is; influenced by distribution of drop sizes, frequency 
of coalescence and breakup of drops, distribution of biomass 
between :fhe continuous phase and the dispersed phase interface, 
mass, transfer between phases, and distribution of substrate in 
the dispersed phase. Rational choice of hydrocarbon-fermentor 
design,. -and, scaie-up has so far been hindered by inadequate 
knowledge of -drop phenomena in these systems. The present 
study .eftempts to thiew some light upon the nature of these 
processes; by proposing a simple model for drop- interactions, 
and eSEplains the drop-size distributions observed in hydrocarbon 
fermentors and similar systems, Further it demonstrates the 
applicability of the model by simulating the drop phenomena 
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and growth in hydrocarbon fermentors on a computer. Coalescence— 
redispersion frequencies needed in the course of simulation have 
been experimentally measured. 

The literature reports various attempts to model the 
growth phenomena in hydrocarbon fermentors. Erickson e^ al . 

(1969 a,b,c,) have proposed models for growth on single hydrocarbon 
drops. These models were used to examine the effects of inter- 
facial area, substrate solubility in the aqueous medium, the 
mass transfer coefficient, and the saturation constant in the 
kinetic equations on the growth behavior. Erickson et al. (1970) 
and Humphrey and Erickson ( 1972)- extended these models to cases 
where swarms of drops exist. These models take into account 
the distribution of drop sizes, the rate of adsorption of cells 
on the drop surface,, the rate of desorption of cells from the 
drop surface, substrate transport between the phases, and growth 
kinetics. The effect of coalescence and redispersion of drops 
was also considered. A discrete uniform distribution and a 
discrete normal distribution obtained from an experimentally 
observed distribution curve ( Chen and Middleman, 19^7) were 
used as drop-size distributions by Erickson et al.(l970) and 
relative importance of each parameter encountered in the system 
was investigated. Drop-size distributions were found to affect 
the batch growth profoundly. It was reported that in the 
presence of a drop-size distribution and little or no drop- 
interaction, the large drops require more time for their 
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substrate to be consumed than the smaller drops, thus increasing 
the time required to complete the fermentation. 

Shah et a^. (1972 a) utilized the model for simulating 
the performance of hydrocarbon fermentors and estimated the 
parameters in the model, ibr parameter estimation, they used 
the experimental data of Prokop et (1971) who studied the 
effects of inoculum size, dispersed phase volume, and substrate 
concentration on the batch growth of Candida lipolytica on a 
model system composed of n-hexadecane dissolved in dewaxed 
gas-oil. It was pointed out that either continuous phase growth 
or growth on small segregated dinps or both, contribute signi- 
ficantly towards the linear growth rate especially when the 
dispersed phase substrate concentration is large. Accurate 
knowledge of drop-size distribution and drop- phenomena was 
emphasized for the sake of better prediction. 

Stamatoudis and Tavlarides (1975) have also highlighted 
the importance of accurate knowledge of drop- size distributions. 
Drop phenomena have been simulated using a Monte- Carlo method 
developed by Zeitlin and Tavlarides (1972). This study which 
utilized the experimental data of Prokop et al . (1971), Moo- 
Young et al. (1971 a,b) , and Wang and Ochoa (1972) for estimation 
of the parameters, concluded that changes of total interfacial 
area with the different drop-size distributions were responsible 
for significant differences in prediction of various size- 
distri'fejution models. Moo-- Young and coworkers (1971 a,b) studied 
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the influence of cell-drop interactions on growth in hydrocarbon 
fermentors. This investigation found that a model which assumed 
that the oil drops were smaller than the yeast cells, gave a 
better fit to the experimental data. The importance of small 
accommodated drops was also advanced by Aiba ejt ( 19 69 a,b), 

1 . 1 Drop-Size Distributions (DSD) in Hydix) carbon Fermentation : 

Recently, a number of experimental studies have reported 
drop-size distribution data measured in hydiocarbon fermentation 
systems. Prokop et ^.(1972), and Prokop and Ludvlk (1973) have 
presented drop-size distributions observed in batch fermentations 
of Candida lipolytica on gas-oil and n-alkanes dissolved in 
dewaxed gas-oil. These observers found that the distributions 
are positively skewed and that a very large number of very small 
drops ( < 3^m) exist. Some distributions, however,, showed a 
bimodal shape at size frequency less than ^io, Distributions of 
drop sizes in samples taken from different positions of large 
vessels did not show any spatial dependence. Dispersed phase 
fraction and inoculum size were found to have a strong influence 
over the drop-size distributions. Distributions taken at various 
times during the course of feimentation were widely different, 
thus emphasizing the influence of biomass concentration on the 
distributions. The physiological state of biomass also seems 
to exercise signifidant effect on the dispersion quality. 
Katinger (1973) reported similar observations from his study of 
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DSDs during femen'ta'bion of gas— oil as well as of n_— liexadecane. 
Blanch and Einsele (1973) measured DSDs during the growth of 
C » tro pical is on pure n— hexadecane "both in circulation— stirring 
as well as in flat- blade turbine systems* Highly skewed distri- 
butions were obtained which were affected by biomass concen- 
tration and the physiological state of the system. Sauter— mean 
diameter in the circulation-stirring systems was not found to 
be well correlated with Weber number. 

Yqshida and Yamada (1971), and Podlech and Borzani (1972) 
measured distributions in similar pure chemical systems. Yoshida 
and Yamada (1971) studied dispersion of kerosene in water in 
gas- bubble columns and in a vessel agitated by a turbine type of 
mixer. Sauter-mean. diameter which did not depend upon sampling 
position, was found to be influenced by power input and dispersed 
phase fraction. Podlech and Borzani (1972) studied distribution 
of diesel oil in distilled water using a radiometric technique. 

An empirical parameter was described to characterize the drop- 
size distributions which were found to be un symmetrical. 

Wang and Ochoa (1972) measured distribution of n-hexadecane 
drop sizes during active growth of Candida intermedia . Both 
the presence of cells and the presence of synthetic detergent, 
were found to influence the specific surface area of the hydro- 
carbon. However, the distribution of droplet diameters was not 
affected by changes in cell concentration above 0,5 grams (dry 
weight) per liter of broth. 
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Scope of the Pr^Rent WhrV ; 

The present work aims at developing a quantitative 
unde..tanai„g of the agglceratl.e and hreafc-„p p^oessea In 
the dispersed ll,uld phase In hydrocarbon fermenters. Thus It 
attempts to resolve a major snag In the theoretical modelling of 
substrate limited growth in these vessels, namely: lack of 

quate knowledge of droj^ phenomena and droi^slze distributions. 

It is proposed to find a suitable model of drop-phenomena and te 
carry out the necessary experimental work to verify the predictions 
of the model, in view of previous experience (reported in 
literature) concerning mathematical intractability of such models, 
it is desired that the model be simple enough so that it can be 
used along with transport and growth processes to predict the 
performance of n-alkane and gas-oil fermentors. This study also 
alms to investigate the effect of various factors, for example, 
dispersed phase mixing, nature of drop interactions, adsorption/ 
sorption kinetics, upon the growth of biomass on lydrocarhona. 

The break-up of subject matter in various chapters is 


as fo How : 

Chapter II cJiscusses the development of techniques for 
the measurement of drop-size distributions in systems with and 
without biomass. Inpirical studies dealing with experimentally 
observed DSDs m hydrocarbon fementors and the correlations, 
have also been presented in this chapter. 
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Coale scence- redispersion model and its implications 
have "been discussed in Chapter III. It includes development 
of the model, its analysis and verification, and the measurement 
of coalescen ce- redispersion frequency. 

Results of computer simulation experiments have heen 
presented in Chapter IV. 

Chapter V summarizes the conclusions and recommendations 
for further research work. 


* * * 
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CHAPTER II 

EXPERIMEtiFTAL ARE MPIBIOAL STUDIES 

This chapter deals with experimental and empirical studies 
of the drop-size distributions In agitated hydrocarbon — water 
systems, with and without biological matter. The first part 
discusses the methods of measuretnen t of distributions of drop 
sizes In liquid— liquid dispersions and their applicability to 
the systems of Interest In this study. In the second part, an 
empirical analysis of drop-size distribution data obtained In 
n-alkane and gas-oll fermenters available from the literature, 

Is described ^d correlations are developed. The purpose of 
this study was to look for certain general trends which may be 
useful for the scale-up of these fermenters, 

2,1 Measurement of Pro id- S ize Distributions (I^Ps) ; 

It Is primarily performed by physical means, e^g. by 
microscopic Inspection, by Coulter counters, by light transmission 
or light scattering techniques, by direct photography, ' and by 
sedimentation analysis. The choice of a particular method 
normally depends upon the dispersed phase system under consi- 
deration and the facilities available. Shah ejfc al .(1972 b) have 
reviewed the various alternatives and their applicability. 

Photography and microscopic Inspection constitute the 
methods most commonly used for this purpose. Shlnnar (19^1) used 
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direct microscopic inspection to measure drop-sisse distributions 
of Shell— wax dispersed in water in a heated vessel. Cooling 
the sample taken out of the system fixed the dispersion, which 
could then be observed under a microscope, and be counted. For 
a system not so easily stabilized, different methods of stabi- 
lization were devised. Madden and McCoy (I 964 ) used an inter- 
facial polycondensation method for this purpose. vanHeuven and 
Hoevenaar (1968) used a concentrated detergent solution to 
stabilize the samples taken from the dispersion. Aiba ^ al . 

(1969 a,b) used direct microscopic observation of samples taken 
from n-alkane fermentors to establish distribution of drop-sizies 
in it. Since the possibility of coalescence could not be 
negated in this process, Yoshida and Yamada (1971) used a flow 
cell adjoining the vessel to photograph ' the dispersions in 
kero sene- water system. Similar arrangement was used by Blanch 
and Einsele (1973) to photograph n-alkane droplets in fermentation. 

Katinger (1970) employed a gelatin embedding method to 
stabilize -samples of dispersions. It consisted of taking a 
small sample directly from the fermentor into a semi-molten 
gelatin solution which was then immediately solidified by dipping 
into ice-water. The drops were prevented from coalescing by the 
high viscosity and the stabilizing effect of solidifying gelatin. 
This method was further improved by Prokop et (1972) to 
measure drop-size distributions in gas- oil fermentation. 
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Direct photography of drops in an agitated system has 
been utilized by many research workers to measure drop-size 
distributions. Kintner £t (1961 ) have presented a review 
of these techniques adopted to bubble and drop research. Brown 
and Pitt (1972) have employed this method for the study of 
kerosene dispersed in water systems. Photographs were taken 
using a photoelectric probe with a synchronized photo-flash. 

Chen and Middleman ( 1967 ), Ward (I 964 ), and Scott al . (1958) 
also followed similar methods for their studies. Sedimentation 
analysis which takes advantage of the specific weight difference 
between, dispersed and continuous phases, was developed by Wang 
and Ochoa (1972) for n-alkane fermentation. A combination of 
sedimentation analysis with radioactive tracer in simulated 
systems was employed by Podlech and Borzani (1972) also. Un- 
fortunately, in biological systems, flocks consisting of hydro- 
carbon droplets, cells and air bubbles, prevent free rise of oil 
droplets due to the buoyancy forces, thus interfering with 
accurate measurement of DSD. On the other hand, microscopic 
methods enable direct measurement of DSD even in the presence 
of these disturbing factors and have, therefore, been selected 
for use in -this investigation. 

Here, the detergen't method ( IM) , the gelatin method ' ( GM) , 
and the encapsulation method (EM) have been compared for their 
applicability to hydrocairbon-water systems with and without 
biomass. During the course of laboratory experimentation in the 
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present investigation, GM and EM were found to be of restricted 
utility , while HI was found applicable to both the systems and 
satisfactory for most of the working conditions. As a result, 
detergent method (IM) was selected for all the subsequent 
experimentation. 

All the experiments were performed in a three-litre 
fermentor. The system was agitated by a Waldhof-type of agitator 
whose details are given in Figure 2.1.1. Dimensions of fermentor 
are given in Table II-1.1, A constant speed motor with facilities 


TABLE II- 1 , 1 s DIMEIT SIGNS OF FEEMESTTOR 


Fermentor volume 

3,000 ml 

Operating volume 

1,500 - 2,000 ml 

Fermentor diameter 

1 28 mm 

Fermentor height • 

255 mm 

Impeller type 

modified Waldhof mixer 

Impeller diameter/fermentor diameter 

0.47 

Height of impeller from bottom 

46 mm 

Impeller height 

78 mm 

Number of baffles 

4 

Baffle height 

55 mm 

Baffle width/ferm enter diameter 

0.1 


to control the speed of shaft rotation from 250 rpm to 1300 rpm 
was utilized. Actual speed of rotation during an experiment was 





caeasured using a stroboscope. In all the experiments, total 
liquid level was initially kept at a level of two litres, and 
therefore an air— liquid interface was always present. Sampling 
was carried out from near the tip of the agitator and the length 
of the sampling tube was kept as small as possible so as. to 
minimize coalescence in the sampling tube# 

Distilled water was used as continuous phase in all the 
experiments. Dispersed phase consisted of either a mixture of 
n-alkane___ or gas-oil. Both the n-alkane mixture and the gas-oil 
were obtained by distillation of a Bomashkin (USSR) crude oil. 
Spectrum of paraffins in n-alkane mixture ranged from Cj to ^2^* 
but the bulk constituted of 0^^ - In gas-oil, the disti- 

llation cut consisted of as much as 20^^ paraffins between to 

^ 24 * 

2.1.1 Detergent Method (IM) ; 

This method was first utilized by vanHeuven and Hoevenaar 
(1968) and is based upon the property of lowering of interfacial 
tension by detergents between the phases. Detergent molecules, 
by virtue of their hydrophilic and hydrophobic groups, form an 
electrostatic film around the drops and thus prevent them from 
coalescing when the stirring has been reduced to a very low 
level or even completely stopped. The dispersion- sample, once 
stabilized, can then be micro photo graphed and counted. 



15 


Some amount of dispersion (10-25 ml) was sampled into 
approximately the same volume of highly concentrated solution 
of Sodium Ibdecyl Sulphate (ShS) which was kept mildly agitated 
in order to prevent coalescence of drops before they are 
completely mixed with the detergent. A drop of the stabilized 
sample was quickly put on a grooved microsco pe— slide and covered 
with a microscope cover-glass, taking care not to retain any 
ai 2 >- bubble below the cover— glass. Once the air— liquid interface 
has thus been eliminated, the distribution could be obtained 
subsequently by counting as it was observed that the distribution 
did not change (indicating that no coalescence of the drops 
occurred) even 4-8 hours after preparing the slide. 

ibr micro photography, the dispersion was placed under a 
Rathenow ROW microscope which was fitted with a camera. With a 
100X magnification, many photographs were taken from different 
portions of the sample on the slide and at different depths of 
focus.* The exposure time for OR WO - RPI 5 film was usually 1/5 
second. The film was developed using OR WO - ED9 fine grain 
developer. Positive prints were prepared on normal Bromofort 
BEH 0 Porte paper for an additional 5X enlargement. Thus the 
total magnification was 500X, and 1 mm on prints corresponded 
to Por measurement of distributions, 1000-1500 drops of 

varying sizes were counted from the prints considering only 
those in sharp focus. 

* 

See Appendix 0 for sample photographs 
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The detergent method could be used for both pure chemical 
as well as fermentation systems. Ibr fermentation systems, it 
was found to be especially good because the microbial cells 
present on the surface of the drops, were displaced from the 
interface in the presence of detergent, and during micro photography 
all the cells were present at a focus different from that for 
drops. It resulted in a more authentic identification of small 
drops. 

2.1 .2 G-elatin- Embedding Technique ( (M) ; 

The stabilization of the dispersions in this method is 
provided by the solidifying-gelatin in which the drops of the 
dispersed phase are embedded. An 8^ gelatin solution is prepared 
and poured into petri-dishes to a 2-3 mm thick layer. Just 
before sampling, these petri-dishes are put in hot water so as to 
melt the gelatin. Sampling is carried out in the petri-dishes 
when the gelatin is just in a semi-molten state. After sampling, 
the petri-dishes are quickly immersed into ethanol-dry ice 
mixture which freezes the gelatin and the dispersion with it. 

The state of gelatin just before the sampling is an important 
factor in deciding the quality of the dispersion sample. If 
the gelatin solution is too fluid, dispersion gets thoroughly 
mixed in it and creates difficulties under the microscope as the 
drops get distributed over a wider thickness thus producing only 
very few dix)ps at any depth of focus. Too fluid a gelatin 
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solution may take a finite time before solidifying and the drops 
have a chance to coalesce during the process. On the other hand, 
if the gelatin gets over-solidified before sampling, the dis- 
persion does not get mixed with gelatin properly and therefore 
the sample does not represent the dispersion in the vessel. Also, 
the petri-dishes must not be kept in ethanol-dry ice mixture for 
too long; otherwise gelatin would crystallize thus obscuring the 
observation of the drops. The rest of the process is same as in 
the detergent stabilization method. 

2.1.5 Encapsulation Method (EM) ; 

The dispersed phase drops present in the dispersion are 
stabilized by coating them with a polymer film. Advantage is 
taken of a fast surface po lycondensation reaction which deposits 
a nylon film at the surface. The reagents used in the study were 
Sebacoyl chloride (soluble in dispersed hydrocarbon phase) and 
Hexamethylene diamine (in the aqueous phase). Sodium hydroxide 
was used as the acid acceptor. Sebacoyl chloride (SC) was added 
to the dispersed phase before mixing it with the aqueous phase. 

The po lycondensation reaction was initiated by adding 50 ml. of 
Hexamethylene diamine (HMDA)/HaOH mixture to the system after the 
mixing was carried out for an appropriate time. After 5-10 seconds, 
the agitation was rapidly brought down from the prevailing rpm 
to 150-250 rpm and a sample was taken. Measurement of drop-size 
distribution was made in a way exactly similar to that in the 
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detergent stabilization method. Reaction mechanisms for these 
polycondensation reactions are discussed in detail by Morgan (196-5) • 
Reaction conditions were adapted to facilitate appropriate inter- 
facial polycondensation keeping in mind the following requirements s 

(a) Hydrolysis of SC dissolved in the hydrocarbon phase 
prior to the reaction is affected by its concentration 
in dispersed phase, speed of rotation and time of 
mixing. The rate of hydrolysis is dependent also upon 
the nature of dispersed phase and that of the aqueous 
medium. Unfortunately, both n-paraffins and the gas- 
oil seem to have a fairly low value of partition 
coefficient (Morgan, 1965) for Sebacoyl chloride and, 
therefore, chloride hydrolysis takes place to a consi- 
derable extent. Tbr hydrocarbon-water systems, the 
hydrolysis was not found to be affected by the presence 
of salts either. This limited the time of mixing 
before reaction to 20-30 minutes, being shorter for 
higher concentrations of SC in dispersed phase. 
Fortunately, this time is sufficient to achieve 
stationary drop-size distributions for small dispersed 
phase fractions. 

(b) Stoichiometry of the reaction of SC and HMDA must be 
found empirically, being a function of dispersed phase 
volume, speed of agitation, and the nature of dispersed 
phase. Especially, the latter affects the partition 
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coefficient of diamine between the two phases 
(Morgan, 1965). The concentration of the acid 
acceptor, NaOH, depends upon the concentration of 
HMDA. 

(c) The conditions of the reaction must be chosen so 

that the polycondensation reaction occurs at the oil- 
water interface and a stable film is formed. It was 
observed that as soon as the agitation is brought 
down to a lower level, sampling of the stabilized 
dispersion into a concentrated detergent solution 
helps to preserve the dispersion even when the polymer 
film formed is weak. 

The reaction of SC with HIM is a very fast one and 
Morgan (19 65) has discussed some of the basic conditions for it. 

The conditions employed in our study were as follows; dispersed 
phase 1-5^ v/v for n-alkanes and 1-1 0?^ v/v for gas-oil; speed 
of rotation 700-1000 rpm; concentration of SC in the dispersed 
phase 0.5 and 1.0/o v/v (higher for large dispersed phase fraction); 
HMDA and UaOH concentrations in water phase 0.35 - 1.0 M and ' - 
0.7 - 2.0 M, respectively. The process of counting the drop- 
size distributions from the stabilized dispersion was the same 
as that for the detergent method. 
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2.1.4 Comparison of Methods ; 

Experimental measurements of DSD in pure chemical and 
in fermentation systems were carried out, Ibr pure chemical 
systems, only IM and IM were employed. Gelatin stabilization 
method could not be used as partial coalescence between dispersed 
phase drops could not be prevented. Samples of dispersions from 
biological systems were stabilized using IM and (M. Eor these 
systems, EM was not employed because of the suspected toxic 
effect of Sebacoyl chloride on the growth of biomass. Various 
distributions have been tabulated in Appendix A. The Sauter- 
mean diameters obtained from different methods during various 
experimental runs are tabulated in Table II-1.2, Reproducibility 
of Sauter-mean diameters was +10,0^, Figures 2,1,2 and 2,1,3 
show the distributions obtained from some typical runs. Figure 
2,1,2 presents a comparison of IM and IM for DSD samples taken 
from a dispersion of gas-oil in water in absence of cells, DSDs 
observed in a gas-oil fermentation experiment are illustrated in 
Figure 2.1,3, It follows from these figures that IM gives results 
fairly similar to those obtained by IM and GM, both .of which can 
be considered as standard. EM and GM have been frequently used 
and reported in the literature (Madden and McCoy, 1964; Mlynek 
and Resniek, 1972; vanHeuven and Hoevenaar, 1968; Katinger, 1970; 
Prokop et , 1972), However, these methods are rather cumber- 
some and of specific utility. On the other hand, IM can be used 
equally well for non-biological as well as biological systems. 
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and is relatively very simple in application. IM becomes very 
handy for measurement of transient state DShs’ v\^here samples 
need to he taken at very frequent intervals and thus results in 
saving a considerable amount of time. Especially in the presence 
of cells, it enables one to obtain better resolutions for small 
droplets, these being clouded by cells in the case of .GM. 

2* 2 Drop-Size Distributions Observed in Hydrocarbon Eermentation ; 

Recently, a number of experimental studies have repo rted 
distribution of drop-sizes in hydrocarbon ferm enters. (P 2 ?okop and 
Ludvik, 1975; Katin ger, 1973; Wang and Ochoa, 1972; Blanch and 
Einsele, 1973)« From the results of these studies, it can be 
concluded that the drop- size distributions in these systems are 
influenced by a variety of factors, such as, dispersed ' phase 
fraction, speed of rotation, biomass concentration, state of 
fermentor broth, nature of dispersed and continuous phases, type 
of agitator, and the geometry of the vessel. The on^ previous 
attempt to theoretically predict these distributions during 
fermentation by Stamatoudis and Tavlarides (1973) using a Monte- 
Oarlo simulation technique involved a prohibitive amount of 
computer time. These authors reported that the exact knowledge 
of DSD is important in predicting the performance of hydrocarbon 
fermentors accurately. Thus a need for an empirical drop-size 
distribution for at least a wide range of expierimental conditions 
is obvious* In the absence of a theoretical model of prediction 
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of DSDs, such an attempt would eventually lead to incorporation 
of real size-distributions into a kinetic model for growth on 
hydrocarbons for extreme cases of complete mixing and complete 
segregation, Theoretical model is, however, still desirable 
in order to take into account the effect of interactions of 
dispersed phase drops which are crucial "^oTf correct p 2 ?ediction 
of growth of biomass (Chapter IV, Section 4.3) in hydrocarbon 
fermentors, 

2*2.1 Normalized Dro-p-Size Distributions ; 

Drop-size distribution data observed in the hydrocarbon 
fermentation as well as from the simulated (i.e. similar pure 
chemical) systems were procured from different sources (Prokop 
,and Ludvik, 1973; Katinger, 1973; Blanch and Einsele, 1973). 

The experimental conditions for which the data were collected 
form a wide spectrum. These studies were performed in the flat- 
blade turbine systems and in the circulation-stirring (i.e. Waldhof 
agitated) systems. The following general trends were obvious; 

(a) The drop-size distributions in all the cases, are 
highly positively skewed. A very large number of 
extremelv fine droplets exist in tbe system. 

(b) Even in the case of gas-oil fermentation, the 
distribution of drop sizes becomes finer with growth 
in the fermentor. 
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(c) The distribution becomes coarser with increase in 
the dispersed phase fraction and finer with the 
speed of agitation. 

( d) The effect of the speed of agitation is more pro- 
nounced in systems with flat- blade turbines than 
in those with Waldhof agitators. 

As has been already observed, the drop-size distributions 
are influenced by a variety of factors which undergo a change 
during the course of fermentation. Figures 2.2.1 to 2.2.3 show 
the drop-size distributions for some experimental runs. In order 
to fit an empirical equation to these distributions, it was 
decided to normalize the DSDs with a suitable parameter thus 
reducing them into a set of normalized distributions. The choice 
of the normalizing parameter is often arbitrary and generally 
a matter of trial and error. Sauter-mean diameter was chosen to 
be the normalizing parameter because it is a mean property and 
can be separately correlated with the experimental parameters. 
Figures 2.2.4 to 2.2.10 present the normalized plots of DSDs 
collected by the different authors. 

Figure 2. 2.1 presents distributions obtained by Prokop 
and Dudvlk (1973) at different times during the course of batch 
fermentation on model gas-oil. Shapes of these distributions 
change in time because of utilization of substrate and build-up 
of cell biomass and extracellular products in Ihe system. As can 
be seen from Figure 2.2.4, these changes can be taken care of by 
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the proposed normalization. During the course of fermentation 
of gas— oil, the changes in dispersed phase fraction are rather 
small and the main contributing factor is the cell density. On 
the other hand, during n-hexadecane fermentation dispersed phase 
fraction also changes appreciably. In figure 2,2*5 normalized 
distributions for this case are also shown. 

furthermore, the distributions are also affected by the 
initial dispersed phase fraction, 0 ^, and the energy input 
(speed of rotation, I), figures 2.2.6 to 2,2.9 present normalized 
distributions for these cases. The data reported in figures 2.2.7 
and 2,2.8 were collected for non- biological (n-hexadecane-- water, 
no biomass) systems using , circulation stirrer and flab-blade 
turbine, respectively (Blanch and Einsele, 1973). The data 
obtained by Katin ger (1973) are shown in figure 2,2.9 and these 
also conform to the previous observations. Also, the data were 
available for DSDs taken by Prokop and Ludvlk (1973) from 
fermenters of different sizes. While the smaller fermenter was 
provided with a single Waldhof-agitator and a draught tube, the 
bigger fermenter was provided with a set of similar stirrers and 
draught tubes, DSDs obtained from these two fermentors showed 
considerable differences, but after nomalization with the 
Sauter-mean diameter they could be reduced to one distribution. 
The normalized plots are presented in figure 2.2.10. The fact 
that similar normalized distributions are obtained in vessels 
which are not only different in sizes but also in their geometry, 
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suggests thai: the Sauter-tnean diameter can serve as- a useful 
criterion for scale-up of fermentors for Single-cell protein 
production. Recently, Sauter-mean diameter has been reported to 
be of use in the scale-up of glutamic acid fermentors from 
n-paraff ins (Hattori et al . , 1974 ). This aspect of applicability 
of Sauter-mean diameter requires further investigation and is 
recommended for further study. 

A polynomial equation of the following nature was found to 
best fit the normalized drop-size distribution data: 

In /fjj(d).d 52 : = % + (d/d 32 ) + ^2 ■ ( 2 . 2 . 1 ) 

i, a/ 

where fjj(d) is the ..number density of drops and d ^2 the Sauter- 

mean diameter, a^^ , b^ and bg are constants and these can be 
different for different systems. Ibr the data of Prokop and 
ludvik ( 1973 ) for model gas-oil feimentation in circulation - 
stirrer (CS) systems, the following constants were obtained at 
5^ level of significance: 

% = 0.58 

b^ =-1.79 -2.05 < b^ <-1.53 

bg = 0.28 0.17 < b2 < 0.39 

Similarly, for several sets of Katinger's data, including those 
in Figure 2.2.9 (covering effects of time, dispersed phase fraction, 
and speed of rotation in non-biological as well as in biological 
systems on model gas— oil .and n-hexadecane) the following constants 


were obtained; 
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= 0.58 

=- 1.86 - 1.96 < b ^ < - 1.76 

^2 = 0.19 < < 0.26 

Blanch and Einsele’s data (Eigures 2.2.5, 2.2.7 and 2.2.8) also 
fit into one equation with the following constants; 

= 0-53 , 

=-1.71 -1.81 < b^ < -1.50 

bg = 0.10 - 0.06 < b2 ^ 0.15 

2.2.2 Correlation of Sau ter— Mean Diameter ; 

Prom the above discussion, it is clear that if the Sauter- 
raean diameter is known, the distribution of drop-sizes in the 
hydrocarbon fermentor can be predicted. Therefore, an attempt 
was made to correlate the Sauter-mean diameter with the experi- 
mental conditions. 

Por the data of Prokop and Ludvik (1973) collected on 
model gas-oil at constant speed of rotation, the following 
equation was found to fit the best; 

d^2 = 0.002 (1 + 3.9 0 ) cm. (2.2.2) 

The mean and standard deviation of the percent error j^defined 

as absolute (<132 prediotea/%2 ejcperlmental ^ ''0°] 

12.01 and 1 1 .‘1 2 respectively. Por data on gas-oil, 

d^2 = 0.002 (1 + 2.7 0 ) cm. . (2.2.3) 

the mean and standard deviation of the percentage error were 
13.54 and 10.47, respectively. 
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The data of Blanch and Einsele (1973) for n-hexadecane 
water system mixed by Waldhof agitator (CS), were found to have 
a very low exponent of Weber number dependence, 

d ^2 - (constant). (2.2*4) 

Katinger's data (1973) including distributions of gas- 
oil or model gas-oil from biological as well as non-bio logical 
CS systems were found to conform to the following correlations 

d ^2 = (constant). ¥e~^*^^ (2.2*5) 

Bor all these calculations, = 5 dynes/cm. was used for 
biological systems, although the value of interfacial tension 
may depend upon the cell density (Prokop et al . 1972j Srivastava 
e t al . , 1970). Por non- bio logical system of gas-oil or model gas- 
oil, = 23*0 dynes/cm. was, used. The effect of 0 could not be 
evaluated properly. 

Attempts to compare these correlations with those reported 
in the literature did not succeed as most of the previous work 
has been carried out in pure chemical systems agitated by flat- 
blade turbines (PBT), Table II-2.1 summarizes the different 
correlations. In most of these cases, the exponent of Weber 
number dependence is -0.6 which is significantly different from 
the ones obtained for the data of Katinger (1973)> and of Blanch 
and EiQsele (1973). In order to check the results, fresh sets 
of experimental data were collected for n-alkanes — water systems 
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mixed, by a WaldJaof agitator. The detergent stabilization method 
had been used fcr the measurement of drop-size distributions. 
Experhnents were conducted at various speeds of rotation and for 
different dispersed phase fractions. These distributions have 
been tabulated in Appendix B. The effect of 0 and speed of 
rotation on the Sauter-mean diameter is demonstrated in Table 
II-2*2# The following correlation :was'‘-foundnv;With_ = 48.1 
dynes/cm, 

d ^2 = 0.0145 (1 + 1.17 0 ) We"°*^^ cm. (2.2.6) 

The mean and the standard deviation of the percent error in 
predictions were 12.0 and 10.36, respectively. Comparisons 
between the experimentally obtained Sauter-mean diameters and 
those calculated using Eq. 2,2,6 are shown in Eigure 2,2.11. It 
can be noted that in this case also, the exponent of Sauter-mean 
diameter dependence is considerably different from the value of 
-0.6 obtained for the EBT systems. A low value of the exponent 
could result from the distribution of power for mixing between 
that consumed for creating turbulence, air-liquid interface in 
case of open vessels, and for circulation in case of CS systems. 
Since a substantial axial flow is required (Brokop and Sobotka, 
1973 ) in hydrocarbon fermentors, a greater attention should be 
devoted to dispersions in aerated systems having a high ratio of 
axial to radial flow. In future investigations, special care 
should be taken in designing a mixer and to the presence of 



43 


0 

p 

B 

D 


O 

to 

o 

O 

O 

P 

* 

• 

• 

• 

• 

• 

0 


<3^ 

C7^ 

VD 

Lr\ 

^1 



CD 

CT^ 

C>“ 

i>“ 

0 


-vl- 

cr- 

CTi 

t>- 

C- 

P 

0 






cvj 


fiicmpa 

b w 

ggBj 

^ PntH 

IZi 

I 5 f 

W pq O 

^ . 1 

]p£j 

mis o 


gMCO 

o W 

!> 

^ P ^ 

o pq P 

P EHiM 

P p w 

p <cj H 

P o:2lg 



1 CO 

O 

CO 

!>- 


ITN 


o- 

T— 



X— 1 

\ * 

• ' 

• 

• 

« 

• 1 

CD 

CM 

o 

OD 

VO 

o I 

1 ^ 

CO 

CM 

CO 

CM 


C-T- 

CM 

CM 

o 1 

0 

cncr> 

r— 

00 

O 1 

•H 

• • 1 

» - 

• 

• 

-P 

KD C- 

1>- 

CJ^ 

LPv 

P 

toto 

CM 

T“ 

T~ 

0 



CM 


LTN 

CM 

m 1 

-r- 

o 

O 

T— 

o 

• 

1 

• 

• 

• 1 

T- 


.r~ 

o 

o 

i 

CM 

CM 

CM 


CM 

ur\ 

IC\ 

,ir\ 

m 

VO 

VO 

O 

O i 

05 

• 

i • 

• 

m 

P 

O 

CM 

VO 

[>- 


to 

CM 


■r~ 



p sn 

°.s 

iq -p 0 
CD d P 
0 -P P 
Ph o •' 
?X3 


O 

O 

o 

o 

o 

O 

o 

o 

o 

o 

VO 

CO 

(J\ 

CM 

LTV 




45 


air-liquid interface with regards to circulation time and scale-up. 

2^3 Effect of Cell Density on the Sau ter— Mean Diameter: 

The study of DSDs in gas-oil fermentors suggests that 
the cell density and the physiological state of the biomass strongly 
influence the distribution of drop sizes in hydrocarbon fermentors 
where all the other’ variables remain almost invariant during the 
fermentation. However, from the existing data no conclusive 
correlation could be derived. Therefore fresh experiments were 
designed to throw some light upon this aspect of the system 
behavior. Experiments were conducted for dispersion of 59^ v/v of 
n-alkanes in aqueous medium in a fermenter fitted with a Waldhof 
agitator. C. lipolytica 4—1 . grown in a chemostat, was added to 
the system in steps. With every step change in the biomass con- 
centration, drop-size distributions were estimated by the deter- 
gent stabilization method after 30 minutes of mixing at 1200 rpm. 

The various distributions have been tabulated in Appendix C, Actual 
biomass concentration at the time of sampling for distributions was 
also estimated by the dry weight method. Sauter-mean diameter has 
been plotted vs. biomass density in Eigure 2.2.12. As was. expected 
there is a considerable drop in Sauter-mean diameter at relatively 
low cell densities, of the oi^er of less than 1 g/l. These results 
corroborate the observations of Wang and Ochoa (19*72) who found 
that the distributions of droplet diameters in presence of 
G. Tvi termed ia cells were very different from those in absence of 
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of cells. An exponential dependence is suggested by the nature 
of the da±a-i Other experiments reveal that this effect may be 
very much influenced by the physiological state of the biomass. 
Por example, there is a tendency for the Sauter-mean diameter, 

increase towards the end of the fermentation in spite of 
the high concentrations of biomass (Prokop and Ludvik, 1973). 
However, this effect could not be studied extensively and is 
suggested for future work. 


* * * 



CHAPTER III 


MODEL FOR DROP PHEHOMEHA 

The perfoimance of hydrocarbon feimentors depends upon 
the extent of dispersion of hydrocarbons in the aqueous medium. 

A good understanding and characterization of the dispersion 
phenomena is essential before rational design of such fermenters 
can be attempted. In this chapter, an attempt towards a reason- 
able quantitative definition of the agglomerative and break-up 
processes of the dispersed phase particles, is described. The 
first part concerns with a coalescence- redispersion (C-R) model 
for drop-size distributions in an agitated vessel and its 
predictions. The second part discusses a method of measurement 
of coalescence- redispersion frequency based upon the C-R model 
and describes the experimental work carried out in this connection. 

3.1 Coalescence- Redispersion (C-R) Model! 

Dispersed phase systems have been a subject of investi- 
gation by chemical engineers for a long time. However, Curl ( 19 ^ 3 ) 
was the first to make an analytical study of mixing of the 
dispersed phase in a liquid-liquid dispersion. !I3iis effort was 
to give a quantitative description to the observations of Rietema 
(1958) on the effect of dispersed phase mixing on conversion 
for various orders of chemical reactions occurring exclusively 
in the dispersed phase drops. Curl’s analysis assumed that all 
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the drops were of identical size and every coalescent interaction 
between two drops was immediately followed by redispersion of the 
parent drop back into two equal- sized drops. The net effect, of 
course, was an exchange of solute between the drops. Since then, 
several attempts have been made to account for drop-mixing. 

Erickson et ( 1970 )have .attempted to answer the obvious objection 
to equal-sized drops by classifying the drop population into size 
categories derived from a discretized version of the observed 
size distribution and allowing for interactions between drops of 
different categories. The authors assumed that two unlike drops 
when coalesced, immediately redispersed into two drops of the 
same sizes as those that coalesced; as before the result of 
coalescence is the exchange of solute material. Obviously such 
an Interaction preserves the size distribution. The assumption 
of two unlike drops coalescing and redispersing into the same pair 
except for the exchange of solute is unreasonable, however. 

Yalentas, Bilous and imundson ( 1966 ) considered both 
coalescence and redispersion as occurring at finite rates and 
predicted drop-size distributions which, however, were not compared 
with the experimental observations. Besides, this framework is 
somewhat difficult mathematically and is not readily applicable 
to the study of rate processes such as heat and mass transfer 
occurring in dispersions. 

More recently, Zeitlin and Tavlarides (1972) have 
simulated the droi>-phenomena in an agitated system by dividing 
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the system into two zones, first the zone of high turbulence 
where breakage of drops predominates, second the zone of circulation 
where essentially coalescence takes place. But no adequate 
information is available about the breakage and the coalescence 
frequencies which appear in the equations. 

The present work aims at developing a coalescence- 
redispersion model capable of predicting experimentally observed 
drop-size distributions in batch stirred vessels, in which lies 
the test of the proposed mechanian of interaction and redispersion. 
The stipulations further require that the model be sufficiently 
simple from a mathematical point of view so that i* is extendable 
to the study of chemical and physical rate processes in dispersions. 
The model features interaction of droplets (distributed according 
to their volumes) with -a certain coalescence frequency and 
immediate redispersion into two droplets of arbitrary sizes subject 
to the constraint of mass conservation. The redistribution is 
assumed to be un if o rm with respect to the parent drop volume in 
the statistical sense. The underlying assumption is that the 
coalescence frequencies are probably of a much smaller order of 
magnitude than the breakage frequencies so that the breakage may 
be regarded as instantaneous relative to coalescence. 

5.1.1 Theory ; 

Let us consider a batch vessel containing two hnmiscible 
liquids, provided with an agitator which disperses one .of the 
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phases as fine droplets in the other phase. The process 
involves an almost instantaneous breakdown of the dispersed 
phase into a large number of droplets leading into a second 
phase ;of droplet - phehomen'a, o_during. whtcfe breakage and coale- 
scence occur at a much slower rate. It is. proposed that soon 
after the second phase is in progress the total number of drops 
remains essentially constant with respect to time. The size 
distribution, however, continues to change as a result of droplet 
interactions which, it is further proposed, results in an exchange 
of mass between the (two) interacting drops. The implication is 
that the drop- in ter act ion is one of exchange of mass rather than 
coalescence and redispersion. 

Let n(v,t)dv represent the number density of droplets 
with volume between v and v+dv at any time t in the vessel. Let 
q(v','^") represent the interaction frequency for the droplet pair 
of sizes V* and v". Given that the interaction has occurred 
between the droplets of sizes v* and v", let the probability of 
a new droplet pair of sizes in the intervals (v’+x, v'+x+dx) and 
(v'’-x, v’'-x-dx) be given by g(v’,v"; x)dx, where x obviously 
varies between -v’ and +v'’. Outside the interval [^-v’ , v''^, g 

vanishes and 
v» 

g(v*,‘V'5 x) dx = 1 (3.1,1) 

Ki.T. KANPUR 
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The population balance equation for n(v,t) is then 
easily derived 


B 


n(v, t) 


1 

2 


q(v’ ,v") 


v’+v"^. V 


g(T' ) + g(v’ 

n(v',t) n(v",t) dv* dv'* 


00 ■ 

- j n(v,t) n(v’,t) dv’ (3.1.2) 

0 

The first term on the right hand side of Equation (3-1.2) 
represents the rate of' formation of drops of volume v due to 
interaction of other drops, while the second is the rate of loss 
of drops of volume v by interaction with drops of other sizes. 


No rate process like mass transfer, heat transfer or chemical 
reaction is assumed to occur in the system. 

At equilibrium. Equation (3.1.2) reduces to 

q(v' ,v'’) |g(v’ ,v''|v-v’ ) + g(v’,v''| v"-v) . 
v’+v"^v ' n*(v’) n*(v") dv' dv" 

CO 

p , ■ * 

' ■ = / q(v,v') n*(v) n*(v') dv’ (3.1.3) 

/ 

j 

o . 

where the terms wi'th asterisk represent equilibrium values. 

At this stage, it is necessary to introduce specific 
assumptions concerning the interaction frequency q(v',v") and 
the function g(v',v"|x) before any prediction of the equilibrium 
drop- size distribution is possible. Let us assume the following 
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type of functional dependence of the interaction frequency upon 
the sizes of interacting dropss 

q(v',v'') = (3.1.4) 

where w^and m are constants. When m is taken to he zero, the 
frequency is constant? a value of ni=l/3 represents the interaction 
frequency ' proportional to the diameters of the two drops; 
m = 2/3 implies a dependence upon the surface areas. 

Let us further assume that the function g(v','v’'l x) is 
a uniform distribution so that 

g(v',v'Sx) = -v’ $ X ^ v" 

= 0 otherwise (3.1 -5) 

With these assumptions, Equation (3.1.3) becomes. 




v’+v" 5 iV 


00 


J 

o 


v’^n*(v) v'‘^n*(v*) dv' 


m. 


( 3 . 1 . 6 ) 


The equilibrium size distribution f*(v) defined by, 


f*(v) = 


n*( v) 


(3.1.7) 


00 

f 


n*(v*) dv' 


may then be ex:pressed in terms of Equation (3.1.6) in the follow- 
ing form by means of a transformation of variables; 
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CD Y 

If- 1 

Y ^ O 


(y-x)“ ^ f*(y-x) f*(x) dx 


QD 


= f*(v) y x“ f*(x) dx (3.1.8) 

o 

If we define (j. *(v) = f*(v), then Equation (3.1.8) 


hecotnes , 
00 y 


^ -y CD 

j i^*(x) dx = jy*(v) f ^*{x) dx 

'TT r\ J J 


(3.1.9) 


which is neadily solved by the tnethod of Laplace transfouns to 
obtain, 


^*(v) = o(.exv ( - ) 

J f 

so that the drop-size distribution emerges as, 

f*(v) = o6v“'^ exp ( - — ^ ) 

The constants c6 and js are gi-ven by 


(3.1.10) 


(3.1.11) 




m / 0 


Ow r- 

Lr-S 

“frr^ J 

(3.1.12) 


[n:3i j Vm 

B 

a 



Eor the case of m — 0, Equation (3.1 .11) predicts an exponential 
distribution for the equilibrium drop- size, 


f*(T) = exp ( - -J* ) 

/^l hi 


(3.1.13) 
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When m = represents the average drop surface and average 

drop diameter for m = l/3* Equation (3*1.11) with m regarded as 
fixed is a one- parameter distribution. 

The .justification for the foms assumed for q(v*-,v'’) and 
g(v*,v "5 x) lies in the ability of Equation (3.1.11) to describe 
experimentally observed drop- size distributions, 

3.1.2 Comparison with Experimental Eata ; 

Drop-size distributions were measured in Waldhof - agitated 
n-alkane - water’ systems using the detergent stabilization 
method. The data obtained for various conditions of speeds of 
rotation and dispersed phase fractions have been tabulated in 
Appendix's. 

Experimental data were also obtained from various sources 
reported in the literature (Shinnar, 1961; Luhning and Sawistowski, 
1971 ? Brown and Pitt, 1972? Prokop and Ludvik, 1973). Data of 
Shinnar, Duhning and Sawistowski, and Brown and Pitt were 
collected in flat- blade turbine systems while those of Prokop and 
Dudvik belonged to the circulation- stirring type. All these data 
were utilized for verifying the predictions of the coalescence- 
redispersion (C-R) model and for comparison of the G-E model with 
other distributions reported in the literature. 

In Pigures 3 . 1.1 to 3.1.6 the drop-size distributions have 
been presented, Pigures 3.1.1 show the data of Shinnar for 59^ 
dispersion of shell-wax in hot water for different speeds of 
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agitation. A continuous line represents Equation 
where the average volume was calculated from the experimental 
data itself, Luhning and Sawistowski * s data in which the 
dispersed phase fraction varied between 0.10 and 0,50, have been 
plotted in Figures 3.1,2. Again, a continuous line represents 
the theoretical predictions using Equation (3,1.13)* Figures 
3 . 1,3 display the data of Brown and Pitt for a range of dispersed 
phase fractions at various agitator speeds. All these data 
show reasonable agreement with the predictions of Equation (3.1,13), 
Figure 3.1.4 represents all' the above data on a single dimension- 
less plot... Ihe scatter observed is perhaps not serious consi- 
dering the diversity of their sources and the errors attendant 
on measurement of drop-size distributions. 

For a more quantitative treatment, a Chi-square test of 
the goodness of fit was perfo rmed . But the test proved too 
stringent ■ being very sensitive even to sporadic appearances of 
one or two *bad data points’. The data were also subjected to 
a Chi-square test using an empirical model of vanHeuven and 
Hoevenaar (19^^) > 


k A 


^A = 


b(k-l)I 


k d 


fc -1 


exp 


where 




1 

( "i “ ^21^ ^3 2 ^ 


and one from Gal-Or ( 1968 ), 

*F( v)= r^P(v’ ) dvV 

o 


k d ' 
b 


(3.1.14) 
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~ d? 

f = 4 — - S exp(-o^d^) (3.1.15) 

L 

T / t 16^^ IS ^2/3 

where oC = ( — ^ j ^ . 

The Chi-square test disqualified these models also for the same 
reason as mentioned earlier. It was, therefore, felt that it 
would be expedient to make a relative assessment of the C-R 
model (Equation 3.1.11 or 3.1.13) with those of vanHeuven and 
Hoevenaar (Equation 3.1.14) and of Gal-Or (Equation 3.1.15)> on 
their applicability to the data at hand. This assessment was 
based upon a simple statistical hypo thesis- testing procedure 
described below. 

Suppose there are two models A and B whose predictions of 
a certain experiment are to be relatively assessed with respect 
to the available experimental observations. Let there be U 

p 

observations available} we calculate the "X values for the qth 

experimental set (obtained by summing for the error sources at 

each point in ,the set) for both models, A and B, represented by 

.2 .2 

X ^ and X 2 respectively. Consider the random variable 

■ i ^ .2 

^ -/-A X.3 • 

The hypothesis, H, to be tested is taken to be 

H ; "I^Model A is either as ^od as or worse than Model B. -J- 
If p is the probability that the random variable X has a negative 
sign, then for H to be true p ^ A/2, Let U_ be Hie number of 
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negative signs of X out of II observations (or sets of experimental 
data). Assuming the observations to be independent of one another, 
the probability of N_ or more negative signs out of denoted by 
P, is given by a Bernoulli distribution as given below, 

( 3 . 1 . 16 ) 


(3.1.17) 

If I is very small, we infer that the observation made 
is very improbable under the hypothesis H, on the basis of which 
we reject the hypothesis and conclude that model A is better than 
model B? the level of significance at which this inference will 
be drawn is given by . 

Calculations were made for different sets of experimental 
observations, The subscripts used were o for Equation (3.1.13), 

V for the expression proposed by vanHeuven and Hoevenaar (Equation 

3.1.14) , and G for that contained in the work of Gal-Or (Equation 

3.1.15) . fable 111-^1.1 displays the results of the calculations. 
It demonstrates the high rating of the exTJonential distribution 
(Equation 3.1.13) over that used by Gal-Or. The model appears 

to be superior to that of vanHeuven and Hoevenaar also. 


P = 


y: (|) p*" (i-p)®-*" 

k=H' 


Since p ^ l/2,we must have 


p ^ y (®) (1)® 


= P 


max 
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TABLE 

III- 1,1; 

OOMPARISOH EETWEEM MODELS 

POR DSL DATA 




'PROM PBT SYSTEMS 


Model 

Model 

To tal 

Humber of 

P 

Percentage level of 

A 

B 

humber 

hegative Signs 

max 

significance at which 



H 

H- 


model A is better 






than model B 

0 

G 

24 

22 

1 .79x10“^ 

0.00179 

0 

V 

24 

16 

0.076 

7.6 


Por tiae drop-size distribution data obtained during the 
present investigation as well as for those of PTOkop and Ludvik 
( 1973 ) j Equation (3<'1.13) was found to give a better fit than the 
models of vanHeuven and Hoevenaar (Equation 3.1.14)» and of Gal-Or 
(Equation 3-1 •15) which among themselves could not be distinguished 
(Table III^r1.2). However, significant improvement in predictions 
was observed when the value of m was different from zero in 
Equation (3.1.11). 

TABLE III-1 .2 s COMPABISOH BETIEEK MODELS EOR LSD liATA 

mm cs SYSTMS 


Model 

A 

Model 

B 

. To tal 
number 
H 

Humber of 
negative signs 

H_ 

P 

max 

Percentage level 
of significance 
at which model 

A is better than 
model B 

Y 

G 

5 

3 

0.50 

50.0 

0 

G 

22 

20 

0. 606x107"^ 

0.0061 . 

1/3 

0 

22 

21 

0.562x10“^ 

0.562x10“^ 
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As a matter of fact, Equation (3.1.11) with m = 0.3 and m = 0.4 
was found to he definitely superior to Equation (3.1.13). Equation 
(3.1.11) with m>0.4 was found to worsen the predictions. However, 
no difference could he ascertained between predictions with m=0.3 
and those with m =0.4. Figures 3.1.5 and 3.1.6 show some of the 
data from Waldhof- agitated systems. Plots of Figure 3.1.5 
correspond to data taken hy Prokop and Ludvfk (1973) during gas- 
oil fermentation. Data in Figures3.1.6 belong to some conditions 
from Appendix B. Phe difference in patterns of these systems 
which have a well pronounced circulation zone, from those mentioned 
earlier which were agitated hy flat- blade turbines must, however, 
be noted. It is this circulation zone which manifests itself 
in terms of higher number frequencies generally obtained experi- 
mentally for large drops than predicted. 

From the above results, it appears that the droplet inter- 
action mechanism suggested here-in has experimental support in 
view of satisfactory agreement with data from diverse dispersed 
systems. The exponential distribution (Equation 3.1.13) describes 
dispersions mixed by turbine agitators. On the other hand, for 
systems having a pronounced circulation zone, the model with 
coalescence frequency dependent upon droplet diameter, produces 
substantially improved results. 
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3 . 2 Measurement of the Interaction greauencv ; 

The present methods for the measurement of coalescence 
frequency can he broadly classified into two categories. The 
first methods involve direct observation of the rate of spread 
of a tracer due to mixing in the dispersed phase (Madden and 
Damerell , 1962? Miller et alcf 1965? Hillestad and Bushton, 1966). 
Por this purpose, some drops having a continuous- phase- insoluble 
dye are introduced into an existing dispersion. Alternatively, 
the existing dispersion may’have’a chemical which on reacting 
with another chemical present in the freshly int 3 ?oduced drops, 
may form a dye. The reaction must, however, be limited to the 
dispersed phase only. In both the cases, it is essential to add 
a dye or a chemical to the dispersed phase, and the contamination 
may result in substantial change in dispersion characteristics. 

In the second method, the speed of agitation is suddenly reduced 
to a lower value and the transient response of an optical probe 
which measures Sauter-mean diameter is noted. The coalescence 
frequency is calculated from the initial slope of the response 
curve. (Howarth, 1967; Blanch and Piechter, 1974). The under- 
lying assumptions are that a step reduction in agitator speed 
results in a step change in the degree of turbulence, and the 
new steady state drop-size distribution is achieved by a coales- 
cence-only process. In practice, however, it is very difficult 
to meet the fiist assumption as some time invariably lapses 



80 


before the turbulence decays to the new value, and the second 
assumption will be valid only for a very specific kind of 
dispersion. Moreover, the use of an optical probe for measuring 
the dispersion characteristics in biological systems is also 
con tiX) versial because of the interfering presence of microbial 
cells and air- bubbles. 

The C-E model proposed in section 3.1 suggests a method 
of measuring the interaction frequency. The method does not 
require addition of any chemical to the system and thus avoids 
any deliberate contamination. Another strong point for the 
method is that the basic model satisfactorily predicts the experi- 
mentally observed DSDs in diverse systems. Either a batch or a 
steady state continuous system can be employed. In the case of 
batch vessels, the method requires measurement of unsteady state 
DSDs at various times from the start of the agitation. According 
to the model predictions, a semi-log plot of “ 2^^) vs 

time should yield a straight line. Prom the slope of this 
straight line, the coalescence or Interaction frequency can be 
calculated. 

3.2.1 Theory ; 

Por the unsteady state case in a stirred batch vessel. 
Equation (5.1 *8) may be written as. 
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00 


n(v’,t) n(v'’,t) 

(v* + v") 


-:g- n(v, t) = 


w„ 


dv’ dT" 


v'+v" >. V 


00 

- Wq n(v,t) !' v’'^ n(vSt) dv’ (3.2*1) 

J 

o 


Sinoe the total number of drops in the vessel is assumed to be 
constant during the secondary phase, 

00 

d , , ff f(v',t) v»® f(v",t) 

f(v,t) = Yi U // — dv’ dv" 

° jJ (v’ + V") 

v'+v" ^v 

CD 

- v“ f(v,t) J v’® f(v’,t) dv’ (3.2.2) 

o 

where n(v,t) =1. f(v,t). 

The taotaent equations may be generated to obtain, 


dt ^0 

= 0 



(3.2.3) 

dt M 

= 0 



(3.2.4) 

— K 

dt ^2 

= w E 

0 1 

r2 2 

rN+1 - 

3 P'm m+2 

(3.2.5) 


Equation (3.2.3) is a reflection of the assumption of constant 
number of droplets in the secondary phase while Equation (3.2.4) 
is a statement of mass conservation. Equation (3.2.5) and those 
of higher moments are an unclosed set when m is not zero. 
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However, for m = 0 (i.e. constant coalescence frequency), 
one obtains, 


d 


P2 


dt 


w N P 

-V (f2- 2^1 ) 


( 3 . 2 . 6 ) 


whicli is easily solved resulting in, 

2 

^"2 ^ ^ 1^1 (constant), exp (- t) (3.2.7) 

so that a semi-log plot of - 2^^^) ys time should yield 

a straight line of slope (-w^ H/3). It would, of course, be 

necessary to test the constancy of and |d^. 

Itefining the coalescence frequency, w, as the number of 

volumes of dispersed phase that undergo mixing in unit time, 

w = w h. 
o 

# 

Hence, Equation (3.2.7) may be written as 

f^2 ~ (constant), exp (- j t) (3,2.8) 

Alternatively, the steady state solution of the continuous 
flow equation may be shomn to be, 


w 


~j~Q ' 


^^/^'m+1 ~P'mhm+2^ 


(3.2.9) 


where P' 2 f second moment of the feed droplet size 

distribution. 


3.2.2 Experimental Work and Results ; 

Transient state ESIte in n-alkane - water system in a 
vessel fitted with a Waldhof- agitator were measured using the 
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detiergenl; stebilizaiiion method. Samples were dravffl from the 
region of intense mixing at various times from the start of the 
agitation. Care was taken so. as not to change the system volume 
appreciably due to repeated sampling. In case large number of 
samples were desired, the experiment was repeated under exactly 
similar conditions. Since it requires the calculation of second 
moments of drop— volume distributions, the experiments were 
carried out under carefully controlled conditions and with high 
precision. In spite of all the care, some samples had still to 
be discarded because some coalescence took place before the 
dispersion could be mixed well with the detergent. Experiments 
were performed for various levels of speed of irtation and dispersed 
phase fraction. The various DSEs measured are tabulated in 
Appendix D. 

In Table III-2. 1 are listed the first moments of DSDs for 
some experimental conditions. After an initial breakdown of 
drops for about 9-12 minutes, the first moment becomes almost 
steady (as predicted by Equation 3.2.4) with a fluctuation of 
+ 30?^ from the mean. It is within the limits of reproducibility 
of our experiments which is obvious if we recall that the Sauter- 
mean diameter in our replicated experiments fluctuates within 
+ 10 . 0 ^. 

Moments of the unsteady state DSDs were calculated and 
2 

( 1^2 ~ ^ P'1 ^ were plotted ys time on semi-log plots, Eigures 
3.2.1 show the plots for various experimental conditions. From 
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TABLE 

III- 2.1: 

FIRST MOMMTS OF DBDP- VOLUME DISTRIBUTIONS 



FROM SOME EXPERIMENTAL RUNS 


time 
(min. ) 


First moment (micron^) 


' 1 5% n-P 
800 rpm 

15^ n-P 

1 200 rpm 

5% n-P 

1200 rpm 

3 

13206 

9102 

8205 

6 

8458 . 

6108 

5312 

9 

5436 

2791 

- 

12 

4324 

- 

2469 

15 

2603 

1150 

- 

18 

1986 

- 

2627 

21 

- 

1 474 

2060 ' 

24 

1671 

- 

- 

27 

2253 

2005 

2080 

30 

2374 

- 

1354 

33 

1648 

1670 

- ■ 


the slopes of these straight lines, values of coalescence 
frequencies (defined as number of dispersed phase volumes that 
undergo mixing in unit time) were calculated using Equation {3«2.8) . 
The results have been plotted in Eigures 3.2.2 and 3.2.3. Figure 
3.2.2 shows that a linear relation exists between w and the 
speed of rotation at constant 0. Similarly, Figure 3.2.3 indicates 
a linear dependence of w upon 0. 
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3.2.3 •Discussions 
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In the present investigation, coalescence frequency has 
been measured in n-alkane - v/ater systems mixed by a Waldhof - 
agitator. Most of the previous work in the measurement of 
coalescence frequencies in dispersed systems has been carried 
out in those agitated by flat-blade turbines (DBT). This has 

been summarized in a review article by Shah et al. ( 1972 b) , 

■% 

These two systems have widely different characteristics. The 
dispersions mixed by Waldhof-agitators are characterized by 
marked zones of high turbulence and of circulation, and are 
therefore referred to as circulation-stirring systems (CS), The 
qualitative observation for DBT systems that w increases with 
mixing velocity and with dispersed phase ratio, holds true in 
CS systems also. Based upon the reported data, it can be 
concluded that for BBT systems, coalescence frequency is propor- 
tional to the powers of N ranging from 1.5 to 3.3, and of 0 
ranging from 0.5 to 1.1. But as is clear from Bigure 3.2.2, w 
in CS systems varies linearly with N. This could be because in 
CS systems' a significant fraction of energy supplied is utilized 
in creation of circulation currents. It has already been reported 
in Chapter II, section 2.2 that in sach systems the Sauter-mean 
diameter is inversely proportional to the 0,1 - 0.3 power of 
the Weber number as against a 0.6 power for BBT systems. Thus 
a rather weak dependence of coalescence frequency upon the speed 
of rotation" is in accordance with other observations on CS systems. 
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However, a careful study of the energy distribution for the 
creation of turbulence, interfacial area, circulation currents, 
and air- bubbles is desired before any generalization can be 
drawn . 

In the present investigation, experiments have been 
perfomed in batch agitated vessels and, therefore, only constant 
coalescence frequency could be calculated* But the assumption 
of the constant coalescence frequency could be done away with 
if the experiments were carried out in continuous flov/ stirred 
tank vessels (CPSTY)* In that case it will be necessary to 
feed the CPSTY with finely dispersed dispersion and Equation 
(3.2.9) can be used to calculate the interaction frequency 
provided the input and output DSDs are known. 


* * * 



CHAPTER IV 


SIMUIATIOH OP HYDROG ARBOR HEMMTATIOH 

In tiiG present chapter, an attempt has been made te 
predict liie perfomance of hydrocarbon fementors using a Monte- 
Carlo simulation method developed by Shah (1974), Both the 
gas-oil and the n- alkane fermentations have been studied. The 
drop- phenomena have been modelled by the coalescence - redis- 
persion mechanism described in Chapter III. Sigaificance of 
dispersed phase mixing and the importance of the nature of 
drop-in ter actions has been danonstrated . The contribution of 
growth in the continuous- phase towards overall growth is also 
investigated, 

4.1 Intro duct ion ; 

A number of attempts have been reported in literature 
to represent substrate- limited growth in hydrocarbon fermentors 
by mathematical models. These models involve various assumptions 
regarding drop-size distributions and their interactions. 

Erickson and coworkers (1969a, b,c) proposed models for growth 
on a single droplet. Extensions to real case fermentations 
were made assuming complete mixing. These models, however, 
present an over-simplified picture of the acinal phenomena 
and were found inadequate to predict the effects of inoculum 
size and dispersed phase fraction as observed experimentally 
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for gas-oil fermentations (Erokop et I971). The model 

for gas-oil fermentation was modified by Erickson _et (1970) 
who took into account the presence of drop-size distribution 
(DSD) and the dynamics of adsorption and desorption of cells 
between the interface and the continuous phase. Dispersed 
phase mixing due to drop interactions was accounted for by a 
mass transfer coefficient term. The interactions pin posed 
were of such a nature so as to preserve the initial drop-size 
distribution which was taken to be the same as that experi- 
mentally observed in pure chemical systems having lean dis- 
persions (Ghen an d Middleman, 19^7). The model is applicable 
only to those systems for which the dispersed phase consists 
of hydrocarbons dissolved in an inert medium. Since it lacks 
the capability In predict the evolution of the drop-size 
distribution, the predicted growth performance is strongly 
influenced by the initial DSD. Moreover, the interactions 
assumed are not realistic. 

More recently, Stamatoudis and Tavlarides (1973) 
have simulated the drop- phenomena in hydrocarbon feimentors 
using a Monte-Carlo method developed by Zeitlin and Tavlarides 
(1972). The biomass growth model is the same as that used by 
Erickson ^ al . (1970) with a slight modification to account 
for reduction in contact surface area between cell aid 
hydrocarbon drop due to rigidity of cell surface. Different 
breakage and coalescence frequencies have been considered in 
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diff'e3?ent parts of tiie farm en tor. The growth patterns for the 
extreme cases of complete mixing and complete segregation of 
dispersed phase drops have heen worked out. However, the cal- 
culations of growth for intermediate cases of mixing between 
drops could not be carried out due to the excessive computation 
times involved. 

In the present investigation, substrate-limited grov^th 
on gas-oil and on n- alkanes has been considered. That is, the 
oxygen concentration in the vessel is maintained above a critical 
level throu^out. It is assumed that the presence of air bubbles 
does not in any v^^ay affect the dr op- phenomena and the gro?>?th. 
Growth of biomass takes place both in the continuous phase as 
well as on the interface between the hydrocarbon and aqueous 
media. Monod's model has been used for representing growth in 
the various phases. In the continuous phase (i.e. aqueous phase), 
the specific growth rate of the biomass is given by, 


S' 

/c .'^max k , + S' 
s' 


( 4 . 1 . 1 ) 


where jj- is the maximum specific growth rate of the microorganism 
/ tn.3,x 

on the given substrate,' k , is the Michaelis-Menten type of 

s 

constant (or saturation constant) , and S' is the dissolved hydro- 
carbon concentration in the aqueous phase. 

The specific growth rate of the cell on n,- alkane droplets 


is asstmed to be, 

V^d max 


( 4 . 1 . 2 ) 


Ibr gas-oil fermentation, the specific growth rate of 
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cells on the interface is assumed to be, 

/^d k + S (4-. 1.3) 

s 

where S is the concentration of the consumable substrate in 
the dispersed phase drops, and is the saturation constant. 

In n-alkane fementation , the dispersed phase itself is 
the consumable substrate. Hence, the dispersed phase fraction 
changes considerably during the course of fermentation which may 
be carried out until all the dispersed phase is consumed. On 
the other hand, only a small fraction of gas-oil constitutes 
the consumable substrate. Therefore, changes of dispersed phase 
fraction in gas-oil fermentation are rather small and can be 
ignored if the concentration of the consumable substrate in 
gas-oil is not very high (Erickson £t , 1969 a). It has been 
assumed in Ihis investigation that the volume of a gas-oil drop 
does not change because of substrate consumption. 

The exchange of biomass between the continuous phase 
and the interface is assuned to be controlled by Langnuir's 
adsorption/ desorption dynamics. Thus, the biomass growth model 
is the same as that used by Humphrey and Erickson (1972). 

The details of the growth equations and the methodology 
of simulation are given in Appendices E and E, respectively. 

4. 2 Simulation of Ir op- Phenomena Without Q-rowth : 

The drop phenomena were represented by the coalescence- 
redispersion model proposed in Chapter III. Any two drops can 
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coalesce with equal parohabil ity ^and ihe coalesced drop redisperses 
iatned iately into two drops of arbitrary sizes. Computer simulation 
experiments were conducted for different initial drop size distri- 
butions (Exponential, Uniform, and Dirac-delta) and for different 
values of 1die coalesceice-redisperslon frequency (CRP). The 
results have been plotted in Eigures 4.2.1 - 4.2.4 and compared 
with the analytical solutions of 'population balance equations* 
(Equation 3.1.13) which have also been presented in these figures. 

Eigure 4.2.1 shows the results for exponen tial initial 
distribution. Cumulative size distributions at zero and at one 
hour of simulation have been plotted. A continuous line shows 
the analytical solution of the population balance equation at 
equilibrium. In the computer simulation experiment, the distri- 
bution of drop sizes does not change with time because the 
analytical solution itself is an exponential distribution. Results 
of un if orm initial distribution have been presented in Figure 
4.2.2, The distribution changes to the exponential form (the 
analytical solution) with the lapse of time. Similarly, also 
for a Eirac- del ta initial distribution, ihe final distribution 
becomes exponential (Figures 4.2,3 and 4.2,4). 

In all the cases, an equilibrium of drop sizes is 
achieved within one hour of the start of the simulation experi- 
ment. Tie actual time for achieving the equilibrium USB, of 
course, depends upon the value of the CRE employed for the 
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experiment. Jbr example, for an experiment mth CRF = 20 Ire”'', 
an exponen tial distribution is established within 30 minutes 
wlereas about double the time is needed for another experiment 
with CRi' = 8 hr”"' (i’igare 4.2.4). 

In order t) determine the size of the system and the 
number of simulations required to produce results within 
acceptable bounds, computer experiments were performed with the 
number of drops equal to 10, 20, 50, and 100. ibr all the 
experiments, 95/^ confidence bands became narrower as the number 
of drops was increased, and also when the number of simulations 
was increased for a fixed number of drops. As an illustration, 
reference is made to Table 17-^2. 1 where percent deviations at 
ijj(v) = 0,5 are recorded for different conditions of the number 
of drops and Ihe number of simulations for Dirac- de 1 ta initial 
DSD. 

TABLE IYT2.1A s 7AEIATI0n OP PER GMT lEVIATION WITH 

H AED E 
os 


Eumber of drops, 

E 

0 

Percait deviation at 
Ej^(v)=0.5 for Eg = 15 

10 

+11.0 

20 

+ 7.0 

50 

i 4.0 

100 

+ 2.4 
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TAjBLE iy-2.1B s YAELATLOS OF PER GMT JEVLATLOE TOTH 

¥ AST) N 
o s 


Number of simulations 


5 

10 

15 

20 


Per cent deivation at 
:^(v)=0.5 forU^ = 20 

± 17-0 
+ 8.0 
+ 7.0 

+ 4.0 


Tbe computer time, on the other hand, increased almost linearly 
with 1he nimber of simulations and exponentially v/ith the number 
of drops. Later, while simulating drop- phenomena v/ith growth, 
it was observed that tte deviations for biomass concentration 
were affected also by the nature of initial drop-si 2 :e distri- 
bution. Por example, for the case of Dirac- delta initial LSD, 
tlffise deviations were within + 1.0 per cent when N and N were 
10 and 15 respectively, whereas for the other two initial USDs 

double the number of drops- v/ere necessary. PorN =10 and Iff =15» 

o s 

simulation of drop- phenomena with growth required one minute of 
computer (IM 7044/1401) time for one hour of simulated fermen- 
tation. Ibr all the subsequent simulations, Iffg was kept at 1 5 
and was either 10 or 20 depending upon initial ISL. 


4,3 Simulation of Drop- Phenomena with Growth ; 

Here the properties of the particles change continuously 
during each IQ and this chaige is governed by the growth equations 
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given in Appendix E, as against simulation without growth where 
properties of the particles do not change except when a random 
event occurs. 

During each IQ, these equations were solved simultaneously 

using Euler's method which was found to he as efficient as a 

Hanmings Eunge-Kutta predictor— corrector method because in 

general, the IQ is very anall. However, the solution of equations 

for substrate concentration in the continuous phase (Equation E.14 

and E.23) required an extremely small stei^-size, Eor instance, 

—5 

for a step size as small as 10 hr., 1iie solution resulted in 
negative substrate concentrations in the continuous phase. Eater, 
an analytical solution of Equation E.14 considering all the other 
variables constant for a very small time interval revealed that 
the substrate concentration ajproaches a steady state value very 
fast, this interval depending upon the value of constant a^. 

(Eor real case fermentations, a^ is of the order of TO.) The 
above problem can be tackled either by using a smaller value of this 
parameter (Ericka>n at ^. , 1970) or by assuming a quasi-steady 
state for the con tinuo us- phase substrate- cone aitration . (Algebraic 
equations E. 16 and E,25 to be used instead of differential 
equations E.14 and E.23, respectively), Using a analler value 
of this parameter is equivalent to increased solubility of 
hydrocarboiis in ttie continuous phase aid results in more overall 
growth rate thai can be expected in actual fermentations. (Growth 
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from 1 g/l to 14 g/l in 10 hours for = 10^ as compared to 
1 g/l to) 10 g/l ?;hen a^ = 10^^,) 

4.3.1- Effect of Segregations in Dispersed Phase ? 

As has heal pointed out by Humphrey and Erickson (1972), 
two types of segregations aie possible. There could be substrate 
concentration differences among drops (possible only witoi gas- 
oil droj®) and cell population differences from drop to drop. 
Differences in substrate concentration arise due to different 
rates of mass transfer from interfaces of drops of various sizes 
and due to different rates of substrate consumption from these 
drops. Substrate con centrati-on differences betv^een drops, in 
turn, cause differences in biomass growth rates on the surface 
of drops, leading to cell population differences. The phenomena 
of coalescence and redispersion of drops, and the adsorption/ 
desorption of microbial cells at Ihe interface tend to reduce the 
degree of segregation. 

In order to study the effect of dispersed phase segre- 
gations, growth perfomances of gas-oil and ^alkane fermentations 
were simulated for the extreme cases of complete mixing and 
complete segregation between drops. These have been presented 
in Figures 4.3.1 - 4.3«4 for two different initial drop-size 
distributions ( un iform and exponen tial ) . In any particul ar 
figure, growth performances for lite extreme cases are signi- 
ficantly different from each other. For example, for the case 
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of uniform initial DSD in gas-oil fementation, as much as 97.5 
per cent fermentation is completed for complete mixing between 
dispersed, phase drops, "while for the case of complete segregation 
only 90 per coat fermentation is over in 24 hours. In terms of 
fermentation time required to convert 95 per cent consumable 
substrate into biomass, a difference of as much as 7 hrs. can 
be expected between the two cases, ilso , these extreme cases 
envelop the growth curves for intermediate degrees of mixing 
discussed below. This establishes that the drop-in teract ions are 
important and cannot as such be neglected while attempting to 
predict the biomass growth performance in these fermentors. 

In order to determine the significance of the nature of 
drop- interactions, the following three types of coalescence-redis- 
persion (C-R) processes have been consideredt 

A. Any two drops can coalesce and toe coalesced drop 
breaks into two drops of the sane sizes as those coalesced. The 
C-R event results in homogenization of substrate concentration 
and the biomass concentration on the interface of the dao-^ps. 

This model of drop-in teract ions can be considered as a modifi- 
cation of Curl's model (1963) and in fact reduces to the original 
model of Curl for the case of Dirac- de 1 ta initial DSD. Since 
the drop- phenomena do not result in creation of nev^ droi^, drop- • 
interactions can be realized in the simulation algorithn discussed 
in Appendix 1, by simply generating the representative sequence 

T.^ 5 U «| 5 T^ > ^2 ^ ^3 ^ U^ 5 • . • • 
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B. The breakage of the coalesced drop is similar to 
that in A but the C-R event does not affect the irdivi*.al biomass 
concentrations on the interface. The substrate mixing is represen- 
ted by a mass- transfer coefficient tem as proposed by Erickson 
M ^.(1970). However, this model of interactions is applicable 
only to gas-oil fermentation. Ibr n- alkane fermentation, this 
model is equivalent to the case of complete segregation of 
hydrocarbon droplets. Realization of this model of interactions 
can be carried out easUy by addition of a term 




in Equation E.17 (Appendix E). 


C. Two drops coalesce and the coalesced drop breaks 
into two drops of ^bitrary sizes. The C-R ev^t results in 
uniform biomass and substrate concentrations in the two drops 
(same as that in A) , This model of interactions has been 
discussed in detail in Chapter lU. 

These models differ fioom each other in the nature of the 
drop- phenomena and the results of the interaction. A schanatic 
representation of the three models is as follows; 
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Computer simulation experiments were cairried out for both 
the gas-oil and the n- alkane fermentations. The value of the 
coalescence-redispersion frequency (CEP) was chosen to be 10-20 hrT''. 
(of the same order of magnitude as that experimentally observed 
in n-alkane -water-¥aldhof- agitator systems. References Chapter 
III), These results have also been presented in Pigures 4.3.1 - 
4.3.4. 

Choice of initial I6D is crucial for models A and B as 
this strongly influences the predicted growth perfoimance. This 
is not the case with model C where an equilibrium BSI) sets- in 
in the early phases of growth during which the growth is small. 

It is only obvious to expect the final growth in actual fermen- 
tations to be unaffected by the initial distribution of drop-sizes, 
since the growth takes place over a period of time incomparably 
larger than that required to establish equilibrium DSD, When 
equilibrium I6D (i.e. exponen tial distribution) is used as 
initial DSD for model A and Model B, all the models assume the 
sane dynamic equilibrium DSD during gas-oil fermentation because 
the sizes of the gas-oil drops have been considered not to change 
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due to substrate consumption. Therefore, any difference in the 
predicted giDv/th perfoimance of gas-oil fermenters will be a 
result of different degrees of substrate and biomass segregations 
in case of ■various interaction models. Interactions proposed 
in model B do not cause mixing of biomass on the interfaces of 
the interacting drops and a high degree of biomass segregation 
is expected. Substrate concentration in drops at different times 
of growth shows the sane trend and, this is why, the predicted 
growth curve for the model B lies nearer to that for the case 
of complete segregation be'fcween drops {Figure 4,3.4). 

Model B represents an attempt to take into account •liie 
effect of dispersed phase mixing by incorporating an additional 
term in the growth equations and, thus, making the problem 
amenable to relatively simple numerical solution. Model A and 
Model C are more realistic realizations of inter-droplet mixing 
and the drop- phenomena, but cannot be easily , subjected to 
numerical solutions. The difference between models A and C lies 
in the redispersion mechanian of coalesced droplets. Thus model' 

0 has the merit of being able to predict drop-size distributions, 
wtereas model A preserves arbitrary drop^size distributions. The 
differaices between predictions of transport and kinetic phenomena 
in dispersions based on models A and C would therefore be expected 
to be due to (i) the role played by the initial size distribution 
and the transient period of evolution ■tso the equilibrium drop-size 



distribution and (ii) the role of the redispersion mechanian 
in distributing the "various component species among the droplets. 
The fiist factor becomes unimportant in dealing with processes 
such as feimentation which are indeed slow in relation to the 
time scale during which equilibrium drop-size distributions are 
attained. The redispersion mechanism is therefore the only 
factor contributing "bD any differences in feimien tation rates 
predicted by models A and C. 

The present simulations indicate little difference be"tween 
models A and C when the equilibrium drop- size distribution is 
substituted into- model A (Figures 4.3.2 and 4.3.4). Thus the 
new redispersion mechanism appears to have contributed little to 
the prediction of f e2?men tation rates in the present context, 
which could not have been asserted a priori . Ho difference 
in fermentation rates can be predicted by the "two models unless 
substantial differences occur in their predictions of the distri- 
bution of substrate and biomass concentrations among "the drop 
population! furthermore large differences in fermentation rates a 
likely "to arise if biomass and substrate concentration ranges 
are such that the growth rate is very sensitive te either one or 
both of them. • 

From the preceding discussion, it should be clear that 
model A cannot always be relied upon to give "the right results. 
Besides, it must be recalled that model C can be used to obtain 
coalescence frequencies from transient batch or steady sta'fce 
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continuous systems of stirred dispersions. Computer time require- 
ments for model A and Model C are comparable when initial DSDs 
are same. In a computer simulation experiment, larger number of 
drops and/or a large number of simulations are required when the 
starting DSD has also to be randomly generated. Drom this point 
of view, model C consumes significantly less computer time as 
all the simulations may be started with a deterministic DSD 
without affecting the grov/th perfomance. Model C, in this 
sense, is a significant improvement over model A. 

The results presented herein establish the importance of 
not only drop- interactions, but also the nature of interactions 
in predicting the perfomance of these fermentors accurately. 
Therefore, in further work, model C has been utilized to simulate 
grov^th performance in intermediate cases of mixing. The fact 
that model G and model A with equilibrium DSD predict similar 
biomass growth perfomance, makes one aware of a new emphasis in 
the areas of interest for further research in hydrocarbon fermen- 
tation, Dor improvements in predictions of growth in these 
fermentors, it will now be desirable to explore the physico- 
chemical and physiological aspects like adsorption and desorption 
of microbial cells on the hydrocarbon - media interface, and 
the mode of uptake of hydrocarbons by the microorganisms. 
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4.5*2 Distribution of DroT>- Sizes ; 

DSDs during the course of computer simulation of gas-oil 
and n-alkane fermentations were also noted down. Since the 
sizes of tbe gas-oil drops have been assumed no t to change due 
to substrate consumption (section 4,l), the distribution of the 
sizes of the gas-oil drops obtained in compiter simulation of 
drop-phenomaaa with biomass growth are the sane as those obtained 
without grov/th. On the other land, the sizes of the n-alkane 
drops decrease with time due to consumption of substrate' and tla 
DSD becOTies finer as the simulation proceeds. However, since Ihe 
time constant for the process of achieving equilibrium DSD due to 
drop- phenomena is much smaller (of the order of fraction of an 
hour) than that for substrate consumption (of the order of hours), 

DSDs in n-alkane fennentations also show a tendency to preserve shape. 
It is evident from Figure 4.3.5 where normalized plots of ISIte 
from a computer simulation experiment of n-alkane fermentation 
have been presented. It may be noted that ihe similar tendencies 
for self preservation of DSD have also been pointed out in 
section 2.2,1 (Chapter II) for experimentally observed DSDs in 
gas-oil and n-alkane fermentations i 

4*3*3 Effect of Growth in Continuous Phase t 

Figures 4.3.6 and 4.3.7 show the effect of mass transfer 
parameter (ag) in simulation of n^-alkane and gas-oil fermentations, 
respectively. The value of this parameter, decides the ex tact of 
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growth 1X1 the continuous phase. A value of 82 = 0.0 when the 
initial aibstrate concentration in the continuous phase is also 
zen, represents the case of no growth in the continuous phase. 

On the other hand, a very high value of a^ implies that the 
continuous phase is almost saturated with the substrate all the 
time. These fi^ires show that the giowth in the continuous phase 
reduces the time required for fermentation considerably. However, 
J'lgure 4.3.6 also shows that the contribution of continuous phase 
growth towards overall growth in n-alkane fermentation for the 
cases of practical, interest (for which value of parameter- 82 is 
of the order of 10“^, corresponding to 1?^ = 9 cm/hr.) is not very 
significant. But, for gas-oil fementations where 1±ie value of 
parameter for = 9 cm/hr.' is of the order of 10”^, very 
significant differences are manifested between the two growth 
curves with and without continuous phase growth. Ibr example, in 
order to consume 80 per cent of the consumable G-substrate from 
gas-oil, as much as 8-9 hr. difference in fernen tation times can 
be expected in presence and in absence of continuous phase growth. 
Thus it seems that for the values of parameter 82 corresponding 
to industrial hydrocarbon fermentations, con tributions of conti- 
nuous phase growth towards overall growth is not as important in 
n-aikane fermentation as in gas-oil fermentation. This fact is 
further supported by the effect of adsorption paraaeter a^ on the 
growth predictions of the two fermentations. These are presented 
in Figures 4.3.8 and 4.3.9. Tbr gas-oil fermai tation, due to 
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higher oontrihution of continuous phase growth, differences 
between the gixjwth predictions for = 1.0 and a^ = 10.0 are 
anall (Figure 4.3.8) as compared to those in n— alkane fermentation 
(Figure 4.3.9)« These results are, however, different from those 
obtained by Erickson et a^, (1970) for gas-oil fermentations. 1 
possible reason f)r this could be Ihe fact that value of parameters 
and a-j^, in present study are 1.0 while those adopted by 
Erickson _et al . were 0.01. The effect of snail values of these 
parameters is to magnify the contributions of both the phases. 
Further studies in this direction are necessary before the contri- 
bution of continuous phase growth can be ascertained for all 
conditions of interest, 

4»3«4 Comparison of Growth Fredictions with Experimental !Data ? 

It was deemed desirable to compare the computer simulated 
growth behavior with that observed in laboratory experiments of 
hydrocarbon fermentation. Choice of writable values of the 
parameters (e.g., maximum specific growth rate, saturation constants, 
adsorption/desorption constants, etc.) that would adequately 
represent the experimental conditions, posed a considerable 
problan, No direct experimental measurement o f these parameters 
in hydrocarbon fermaitation systems has yet been , reported; even 
solubility of long-chain hydrocarbons in water is a matter of 
controversy (Mcluliffe, 1969? Yoshida and Yamane, 1971 j Bell , 1973; 
Goma, ^ al, , 1973; Baj pai and Pro kop, 1974) , 
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Shah et ( 1972j^_)estimated the paraneters hy minimizing 
the sum of error squares between the predicted and experimentally 
observed growth data in gas-oil fermentors. Stanatoudis and 
Tavlarides (1973) also adjusted the parameter values so as to 
fit the predicted perfomance into experimaital data. Even though 
tie experimental data utilized by the above groups of worke^rs are 
the same, considerably different values of the parameters were 
obtained, Ibr exanple, maximum specific growth rate, n reported 

jT m S.X 

in the first case is 0.75 as compared to a value of 0,25 in the 
secon:d case; (as well as k^) values are alsD different by two 
orders of magnitude. After a considerable del iberatiDn, the 
following parameter values were chosen; 

Aiax = 

ka = 0.166 liter^p /(hr., gi Icellp^,) 

k^ = 0,61 hr””* 

kg = 16.6 g/literpp 

kg,= 3.07 X 10”^ g/liter^p 
Seq* = 4,9 X 10”*^ g/literQp 
n^^ = 8,16 X 10“^ g cell dry wt./cm^ DP 

These terms have been explained in Appendix E. It may be noted 
that these values are same as liio se obtained by Shah et (l972fXi^ 
for one of their models, with a minor adjustment in the value of 
/*max* ^3r)v?th data obtained by computer simulation are plotted 
in Pigsres 4,3.10 and 4,3*11 • Average volume of gas-oil drops 









Mode-; gas-oil form< 2 nt< 

Obsorvod growth; dato,; 
o^xperifT'ent'# 39 ■ ' ^ 

( PTok c p gt Qi ^ 137 1 ) : ' 

"heor 2 liC''diy predicted; i 
curve using 'modei C’; 


5o=0-06-4h 







127 


at the start of simulation was calculated using DBI) data observed 
by Piokop and coworkers (1971, 1972) for the experimental conditions. 
There is a reasonable agreement betvi'een predicted aid experimental 
results. 

However, a word of caution must be introduced regarding 
the values of parameters chosen. This set is in no way unique 
and it may be possible to find a better set of values. In fact, 
parameter values obtained in this manner normally show a sigii- 
f leant bias towards the model chosen and may not represent the 
actual values for the system, v;hich must be determined by indepen- 
dent measurements. Extensive studies, both theoretical and 
experimental, are needed to achieve this and are recommended 
for future work. 

In the above pages, an investigation into the effects of 
various factors upon the growth in gas-oil and in n-alkane 
fermentors has been reported. The coalescence-redispersion model 
(discussed in Chapter III) is a realistic representation of drop- 
phenomena, which is necessary to predict the performaice of these 
vessels accurately. This model is expected to perform better 
than those previously proposed because of the experimental support 
various components of the model enjoy. 


* * * 



CHAPTER Y 


CONCLUSIONS MD BEOOMMIHPATIOMS 

A coaLe scGiice— redispersion (C— R) model has been proposed 
for the drop- in ter act ions in stirred liquid-liquid dispersions. 

It is a simplifitjd representation of the complex drop- phenomena 
taking place in such systems and yields an analytical solution 
for the case of equilibrium drop- size distribution (LSI)) in batch 
vessels. Its capability to represent a large number of experi- 
mentally observed LSL data obtained from diverse sources (including 
those taken during the present investigation), gives credence to 
the mode of interactions proposed in this study. The model predicts 
a single parameter expression for equilibrium DSL in stirred 
vessels. This parameter can be correlated to the experimental 
variables and, then, the DSD corresponding to the experimental 
conditions can be theoretically predicted. 

Experimental measurement of DSDs in n-alkane-water-Waldhof 
agitated systans, have been carried out using a detergent stabi- 
lization method (IM) developed in this study. This method was 
found to be especially suitable for studies in biological systems 
and is relatively simple. - 

DSIte observed in fementation systems, have been f)und to 
be self preserving in nature. Sauteii-mean diameter, ihe normali- 
zing parameter, takes into account the eifect of dispersed phase 
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fraction, speed of rotation and state of feroentor broth. It 
When correlated with exjertaental oonditlong, revealed that the 

circulation stirring systms behave in a way different from the 

flat- blade furbine systems -Po^+ ^ 

j - • ills fact was also coric bo rated by 

the correlations of 0-E freciuency with tte experimental conditions. 

A new method for the measurement of coalescence- redlsperslon 

frequ®cy has been proposed, based upon the 0-R model. Constant 

scene, frequency has been measured by carrying out studies in 

batch stirred vessels and has been correlated to the experimental 

conditions, ifcwever, sine depadent C-E frequency can also be 

measured by using a continuous flow stirred tank vessel ( CESIV) , 

and is recommended for further work. ITse of OPSW will make the 

measurement quite easy as only two distribution - sampl es, those of 

incoming and outgoing streams In steady state operation, need to 

he taken as compared to many distributions which must be measured 

while working with batch vessels. Moreover, replicate saiples 

can also be takai with as much ease. 


in attempt to simulate the perfomance of hydrocarbon 
fementors using the C-R model along with a model for growth, of 
biomass in the system, was made using a Monte- Carlo simulation 
procedure. The importance of the extent and the nature of drop- 
interactions, and contribution of continuous-ptese giowth towards 
overall growth have been investigated. The simulation confirms 
the self preserving nature of DSDs, observed experimentally , 



in hydrocarbon fermentors. The computer time required for tie 
simulations was also within reasonable limits. This simulation 
procedure can be easily extended to determine tile optimum strategy 
for the operation of batch hydrocarbon fermentors. It would, 
however, require additional information regarding the cost of 
raw materials and the proteins, operating costs, aid "file recovery 
costs of sin^e cell proteins and the residual n-alkanes. 

Aithou^ this Investigation throws some light upon the 
nature of drop-size distributions in hydrocarbon fermoitors, 
further studies ( bo th theorBti.o41 as well as experimental) are 
desirable ti determine the kinetic importance of various sizes 
of drops. It has already been advocated in literature that the 
submicro SCO pic drops which could be aggregates of dispersed phase 
molecules play a crucial role in deciding the mode of uptake of 
hydrocarbons by the microbial cells. This brings us also to the 
question of solubility of hydrocarbons in aqueous phase about which 
many conflicting data have been reported. Theoretical studies 
supported by experimental work in this field, are needed to 
explain the diversity of aa ch experimental findings. 

lastly, the knowledge regarding distribution of biomass 
between the interface and the continuous phase, is a major 
missing link in modelling o f hydro carbon feimentors. So far, 
it has been tackled only by hypothesizing the phenomena. A 
fhndanental inrestigation in to the nature of interactions between 
cell surface and lydro carbons, and cells aad aqueous phase is 
desirable. 


* * * 
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appmdij a 

PBOP-SIZE DISTHEBUTIOES MEASCJREP USUff 2EEERGEIT 
MEHiOD (IM), GELAOIE MEEBOP (CM), AHD MCiiPSULATIOF 
METHOD (M) OP STABILIZATEOII OP USPERSTOlf SMPDES 
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TABLE A-1 s BISBERSIOL SAMPLED lURIffG CQHTINnogs GBL'PIVATIOl 

OE G^LIPOLYTICA OH gAS-QIL 


Dispersed phase fraction 
Type of agitator 
Speed of rotation 
Biomass concentration 


Experiment No.1 


0.10 
l^ldhof 
1000 rpm 
no t noted 


Sl- 

Bo. 


Size Category 

( 


Number of Drops Observed 

m (M 


(immediately (1 hr after 
after sampling) sampling) 


1 

2-4 

661 

943 

479 

2 

4—6 

218 

459 

197 

3 

6—8 

104 

235 

68 

4 

Sr12 

89 

201 

94 

5 

12rl6 

60 

86 

84 

6 

16-20 

44 

59 

70 

7 

20-24 

41 

48 

62 

8 

24-32 

68 

81 

93 

9 

32-40 

25 

35 

44 

10 

40-48 

17 

22 

38 


Total 


1327 


2169 


1247 
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a SIEBSIOH 0? P-ALMES M TIT RTTTt.t;t. 


Dispersed phase fraction 

Tjpe of agitator 

Speed of rotation 

Concentration of Sabocojl 
chloride in dispersed 
phas e 


0.04 

Six-hladed turbine 
700 rpn 

0.5?^ v/v 


Size Category 
( 


jrumber of Drons ObffOTnmd 


8-10 

10-12 

12-14 

14-16 

16-20 

20-24 

24-28 

28-32 

32-40 


To tal 


1675 


1751 
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TiBLE A m SPERSIOl OF GAS-O IL II DISTILLED WiiTOP 


Dispersed phase fraction 

Type of agitator 

Speed of rotation 

Concentration of Sabacoyl 
chloride in dispersed 
phase 


Experiment Wn. e; 

0.08 

six-bladed turbine 
1 000 rpm 


1?^ v/v 


SI, 

No. 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 


Size Category 

( /im) 


Humber of Drops Ob.c!(=‘ 3 ;»ve( 3 , 


0-2 

2-4 

4—6 

6-8 

8-10 

10-12 

12-16 

16-20 

20-24 

24-32 

32-48 

48-64 


Total 


IM 


. 1030 
275 
141 
62 

38 
26 

39 
12 
10 

5 

8 

1 


M 


785 

272 

174 

70 

16 

8 

11 

8 

3 

1 

1 

1 


1647 


1350 




gSPERSM SMPLED lURIHG BATCH CULTIVATION 
01 C, LIIOLYTICA ON GAS-OTT. 


Ex-periment No. 7 
Dispersed phase fraction ; 0.05 

Type of agitator : Waldhof 

Speed of 2x>tation s 1000 rpm 

Biomass concentration ; i g/i dry wei^t 


Size Category 

m) 

4—6 

6^8 

8-10 

10-12 

12-16 

16—20 

20-24 

24-28 

28-32 

32-40 

40r48 

48-56 

56-64 

64-72 

72-80 


Number 

_pf Drops Observed 

IM 

m. 

190 

273 

66 

125 

40 

107 

30 

66 

41 

69 

23 

56 

22 

32 

13 

30 

6 

15 

21 

29 

13 

16 

5 

14 

1 

7 

2 

3 

0 

3 


!&> tal 


473 


845 
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PISgEBSIOIT OP SAS-OIL m T1TSTIT,L1!!T1 WiTTO 

Experiment Nn«, g and in 

Dispersed phase fraction : 0.10 

Type of agitator : modified Waldhof 

Speed of rotation s iooo rpm 

Concentration of Sabacoyl 
chloride in dispersed 


SI. 

No. 

Size Category 
m) 

Number of Drops Observed 

m 

(expt.No.9) 

EM 

(e 2 :pt.No. 9 ) 

IM 

(exp. No. 10) 

1 

0-2 

673 

840 

1310 

2 

2-4 

347 

300 

429 

3 

4— S 

204 

133 

136 

4 

6-8 

109 

102 

86 

5 

8-10 

64 

67 

37 

6 

10-12 

40 

26 

22 

7 

12-16 

83 

58 

32 

8 

16-20 

51 

48 

12 

9 

20-24 

38 

20 

7 

10 

24-28 

35 

13 

8 

11 

28-32 

21 

10 

6 

12 

32-40 

27 

9 

7 

13 

40-48 

11 

4 

7 

14 

48-56 

8 

0 

2 

15 

56-64 

1 

0 

0 

16 

64-72 

0 

0 

0 

17 

72-80 

1 

0 

0 


Total 

1713 

1630 

2101 
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^S.?BRSIOU Qgn- AIMES IK DIFiT-rT,T,T^.T^ garpm. 


Type of agitator 

Speed of rotation 

Concentration of Sabacoyl 
chloride in dispersed 
phase 


Experiment Ho. 1 1 


Dispersed phase fraction ; 0.08 


Waldhof (modified) 
1000 rpm 

1?^ t/v 


(Samples collected ^ter 20 minutes from the start of 

agitationT 


Size Category 

( /im) 


Humber of Drops Observed 


2-4 

4-6 

6-8 

8-10 

10-12 

12-16 

16-20 

20-24 

24-28 

28-32 

32r40 

40-48 

48-56 

56-64 

64-72 

72-80 

80-90 


Total 


1142 


947 
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-^-7 = BlgEERSroH OP d.JLEMES hi TIISTTT.T.CT w.inri, 


Experiment Wa. 16 


Dispersed phase fraction : 0.0^5 

Type of agitator j modified Waldhof 

Speed of rotation : lOOO rpm 

Concentration of Sabacoyl 
chloride in dispersed 
phase . 

(Samples collected after 5^ minutes from the start of agitation) 

SI 

• Size Category 

Number of 

Drops Observed 

No 

• ( /I m) 

IM 


M 

1 

0-2 

1263 


1438 

2 

2-4 

. 233 


213 

3 

4-6 

' 77 


73 

4 

6-8 

40 


34 

5 

8-10 

41 


24 

6 

10-12 

20 


21 

7 

12-14 

11 


12 

8 

14-16 

16 


12 

9 

16-18 

7 


5 

10 

18-20 

10 


4 

11 

20-24 

5 


13 

12 

24-28 

6 


8 

13 

28-32 

8 


2 

U 

32-36 

1 


2 

15 

36-40 

0 


1 

16 

40-48 

2 


1 

17 

48-56 

1 


0 


To tal 

1741 


1863 
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APPMIOIZ B 

EPHDCT OP DISriRSED PHASE PRACPIOH AMD SPEED OP ROTAPTOH 
OP MIXER UPON THE DISTRIHJTION OPn— ALKJJE DROPS IN DCS- 
TILLED WATER IN A WALDHOP- AGITATED PBIMMTOR. 



MM Brlj mSPERSION OF 2io n- ALKANES IN DISTILLED WATER USING M0DIPIE3) V/ilDHOP MIXER 
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APPMDIX J) 

DISTHTHITOT OP n-iLME lEOPS DISPERSED IH nESUDLED 
WAOIER BY A WADDHOP IGITAODR AT DEPPEEMT TEMiS iE(M 

the start op MIHBG 



TABLE D-1 ; DISPERSION OF n~ ALKANES IN EESTILLED WATER 

Dispersed phase fraction: 0.05 
Type of agitator: modified Waldhof 
Speed of rotation: 800 rpm 
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TABLE D-5 ! DISPERSION OE n- ALKMES IN DISTILLED WATER 


T 


a 

Pi 

p 


m <N 
-H o 

O -H *Td O 

^ ^ r -4 O 

• O cs! o 

O' s ^ 


o 

•H 

43 . 

o 

05 

U 


05 


^ 05 
© 

m ch 
P o 
© 

ft © 
© ft 

*H i>:j 

ft ^ 


SO 

o 

•Hi 

-P 

05 

-p 


ft 

o 

rd 

© 

© 

ft 

m 


u 

©I 

© 

ft 

lO' 

m. 

Pi 

I 

M 

Oi 

p 

© 

[ft 

Pi 

Ift 


ft 

fft OD 
K >00 


tP 
K\ 


O 

o co^ 

K\a 3 


CM M 


fti 


pq 

(CM O 

•H 


0 

01 

VD 


c o 

M Jsj 

a 

a ^ 
|ft ^ 

:e^ m 


§) 

©. 

•p S' 

05 

o 


© 

tS 3 

•H 

CD 


rH, 0| 


i 


00 ^ C3^ CO ^ O- CM' K\ o CTv ^ T- ^ | 
i>'’r-cM ir\ moD 00 VO 
cm: T~ T- 


'H- tn K^VD CT^ t>-tOiO tc\ 00 r- in tA I 

tnino tncno-o-vo ^-=;Mn cncM 

V*** 'V~“ 


CO^lC^OOI^“CMOOOVOO^-OT- mtn CM I 

sncoocMC>~<T\£>-inEninK%cM^ 

CM CM CM T-“ T- 


tnCM ^I>-CM T- t^CM £>-CM 00 T- CTvVO I 

^ 001^^00 t<^CMo^£>“lC% 00 l>-CM’r~ 


tn in VO "^cnT-'^tnoT-'£>-co ir— voc^cm 
cn c>-vo T- T- cnco cn covo go vo '^tncM «r- 


■^c^o-c^o-vovot--^ cr>vo mo o im ^ -c- 

mtn'H'mmcMT-CMCMCMVD i>-m^T-- 


mmvoovovo CMr-o-mvo mcMCMCMcMin 
oO"^^mmcMCMCMCMtnvot>-covovDvoT“ 


CM ^VO 00 
O CM 


O CM ^ VO GO O CD CM VO O 00 VO 
T-T-T~T-T-'CMCMCMtnm*^'«M-in 
II i I 1 f i I t I I I I 
CO o VO 00 o 00 ^ VO O CO 

T-T-T-T-T-CMCMCMtncn-^^ 


T“ CM tA'^ mKo c-~oo cno 't- cm tn ^ invo {>- 



VO 

CM 

VO 


in 

O 

cn 


o 


CM 


o 

o 

tn 


CM 


cr\ 

in 


in 

S 


rH 

05 


s 


156 




TABLE D-6t DISPERSION OF n~ALKMES IN DISTILLED WATER 
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TABLE D-T i DISPERSION OF n~ ALKANES IN DISTILLED WATER 

Dispersed phase fraction j 0.15 _ 

Type of agitator 5 modified Waldhof 

Speed of rotation s 800 rpn 
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APPMECX E 

POPULATION MD GROWTH EQUATIONS 

Population Balance Equation : 

A droplet present in the n-alkane fermentation system 
is characterized by its size and the amount of biomass it has 
on its surface at any time. Hence, the population density of 
drops can be considered a bivariate ftinction, the variables being 
V (the volume) and n (the density of biomass on its surface, 
defined as the amount of biomass per unit area). 

The drop- phenomena is considered to be represented by the 
CO ale sconce- redispe 3^ion (C-E) model discussed in Chapter III. 
Briefly, any tv/o drops present in the system are capable of 
coalescing v?ith a constant frequency, w. The coalescence is 
immediately followed by redispersion into two daughter droplets 
whose sizes are uniformly distributed in statistical sense. The 
C-R process is assumed to be instantaneous and as a result of 
this, surface biomass densities on the two drops become uniform. 
It is also assumed that the drop - phenomena do not cause any 
exchange of biomass between the two phases and that these do not 
interfere with mass transfer process. The population balance 
equation for the coalescence-redispersion model can be easily 


written as 
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f(v,n,ttf + ~ f f(v,n,t) 


dt 


+ ^ |f(v,n,t) 


drt 

dt 



n 


1 max 


dv^ dr^ 




J dn^w gCv^ ,T2 ;v) f(v^,n^,t) f(v2,n2»'t) 


a> n 

CO 

- J dY^j dn’ 


w f(T,n, t) fCV-jiiSt) 


(E.1) 


o o 


The first tem on the right hand side corresponds to the rate of 
gooeration of drops due to the 'C-R events, and tte second term 
represents the rate of loss of drops under consideration, due to 
coalescence with other drops. f(v,n,t) is the population density 
of drops, and f(v,n,t) dv dn represents the fractional number of 
droplets at any time t having volumes between v and v+dv, and 
biomass density between n- and n+dn, n2 is the surface biomass 
density on a drop of ■volume V2 which after coalescitig with a drop 
of 'Volume v^ and airface biomass density n^ , and redispersion will 
produce two drops having surface biomass density n. By biomass 
balance, n2 can be expressed as 


f^^2/3 V 2/5 

“2 - “ V2^ + - v^^ 


V. 2/3 




(E. 2) 


max maximum surface biomass density on the drop of 

vo lume v^ that is capable of producing after a coalescence- 
redispersion event with a drop of volume V2, drops of aarface 
biomass density n; it can be easily shown to be, 
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n 


1 msix 


= n\ {^ ) 
1 / 


2/5 


+ (: 


+ 


1 - ^ ) 
"^1 


2/3 


(E.3) 


Population balance equatioiis for the ease of gas-oil systems can 
be written in a similar way. The population density of drops 
in this c^evdll be a trivariate function of drop volume v, 
surface biomass density n, and concentration of conajmable substrate 
S in the drop. 


Growth Equations for n-JO.kane Pemen tation : 


Ibr a single droplet 

n.A - ijj (Sg^,-S').A (E.4) 

ft 

where A is the interfacial area of a drop of volume v, n is 

CO 

tlB maximum surface biomass concentration. Tie quantities with 
primes correspond to the continuous phase. 

Material balance for biomass and aibstrate in the continuous 
phase leads to: 


CO n 


s 


o o 


(E.6) 


f(v, n,t) dn dv 

■ft =Aiax VTTp- - J I ta ^ 

® 0 0 


fCvjn^t) dn dv 


(E.7) 
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If Nq is total number of drcpopresent in tlie system 
and a subscript i denotes each drop, discretized forms of Equation 
(E.4) to (E.7) can be written as 


It ( f'sn) = \ 

-y- i = 1 j 2 y * • . 


dt 


“l) =A,a3j Vl ('’cD-‘’l>°'V'"d “U-l 


Y i = 1,2,...,!^ 


(E.8) 


(E.9) 


4 /'tr . N Z' max 


N 


at^'^o s') =- 


0 


^ +'S' ^ \i ^ 


s 


i=1 


II 


(E.10) 


•|:r( V n » ) = n n'V -y I 

dV c Ama* k„, + s c 2_ ] 

® i=1 ^ 


k n’ (n -n . ) - k 


CO " 1 ' "d^ij "S. 

(S.11) 


On non-dimensionalizing these equations we get, 


if^i = - ^i a^d-x*) x^2/3 

-V-i=1,2,. .,11^ 

d 2 

dj^i ~ ^i y' ~ ^ W” 

-V- i=1,2,. ..,H 


dj 


’ f — 




1 

^ aj y + aga^ (1-x’) 


i=1 


F 


o 


X 


I- y* = 

ajc’ 


)_ X. 

i=1 


.2/3 


(E. 12) 


(E.13) 

(E.14) 


a4(i-yi)y ’-a^yi( e.i 5) 
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where 


1 V n ’^~S ^ y = Y~i » /= t ;a 

c OD , ^ fo Am 


Kt 

CO , ^rn 

1 '“ ■ '"’•t Z’? ’ * ^■p — — — 

Y f ^ n 

« c Amax 


^eq ' 


u y ”'/5 P 
/max '^c ® 


a 


^tr = ^ 

^eq • 


, c6 = ( 36 tt) 


Assuming a quasi-steady state for substrate concentration in the 


continuous nhase 


-B + + 4aj^' 


(E,16) 


where B = _ 1+ a, , + 


^2 r 


Unless stated otiierwise, values of 1hese parameters used in 
simulations^ were: 

Ag = 0.85 g/caP 

substrate 

= average volume of drops at the start of simulation 


= 1 X 10~® cm^ 


0 = 0.10 


“o “ initial surface biomass density = 0.0 gm/cm' 


^co ~ 2 X 10 gm/cm^ 


Aax = °-'75 


kg , = 2 X 10‘ 


gn substrate/L, 
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“■ 0»015 l-Qjp / gn 


= 0.75 ]ar“ 


"o = ^eq» = 2x10" 


@n substrate/L, 


Y _ 0.85 gpi cell^^/gn substrate 
= 9 cm/hr 
GRP = 10 hr""' 

= 10 or 20 depending upon vshether initial BSD is Dirac 
delta or uniform (and exponential), 

N3=15 

Growth Equations for Gas-Oil Pemen tation ; 

Equations corresponding to (E.8) to (E.ll) are 


^ c! ^ ^ max , , , 

at = - — ip-g, \ h - \ (Se,, -S') \ 


^i=i,2,...,ir^ 


(E.17) 


dt ^i^ /Wax k_+S- °i\ 


V-i=1,2,.. 


(E.18) 


( Y 
dt 


ju ^ 

(v et 1 _ /Wax §* ^ , , , , 

^ c' S ) - Y kg.+S’ "" ^ Y. ^i^^eq! 


(E.19) 


« =Aa.E;7?5' - Z - Vl} 

(E.20) 

where is the continuous phase substrate concentration in 

equilibrium with drop i, and is expressed as 


Son * 
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On non- dims 


sionalizlng these equations, we get 


= - ^1 ^i‘ 


- a^Cl-y’) Z.-V3 

» * • • ( E,21) 


H - ^3 jj 


-¥-i=1,,. ..,51^ .(E.22) 

d f ^ 

^yX y» + a2a^ ^ (e.23) 

^ a^.+3f' ^1 L ^i I a^y'(l-y^)-a^j^l (E. 24) 

s ^ i 


where 


y '^Aaax’ ^i - 


n^ V. 

’ ^i “ rT" ’ ^i ~ Y~ ’ 


X* = — , y, ^ __nl. 

Seq* ’ ^ Y S 


Q — m • t - 

^ 7^3 S„ ’ 

/max c /-UK 


a„ = 
4 


Y S 

8 > O 


aj- = -2. , a = ^ a, - 

5 s^eq*’ ^ ^0 ’ "i;;. 


Assuming a quasi-steady state fcr substrate concentration in 
the continuous phase 

-BH.i/9 + 4a,^,.q 

2 ^ (E.25) 

Where B = -,o+ a^, . — ^ ; a = "" 

^2 y ^ X° 

i=1 i=1 
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Unless stated otherwise, values of these parameters used in 
simulations were: 


= 0.85 gn 


/coi' 


substrate 

®°substrate/c®l)p 
1 X 10~® cm^ 


0 = 0.10 


= 0.0 ga cellj^/om^ 

^CD = X 10“^ gn cell^^cm^ 

Aax = 0-75 hr-'' 

^s' = 2 X 10 ^ gn su bstrate/L^p 

= 0.17 gn su bs trat e/ccipp 

\ = 0.013 Loj/gn cellp^hr. 

= 0.75 hr"'' 

I , _6 

^ “ ®eq' = 2 X 10 gn subetrate/Lj^p 

Y = 0.85 gn cell™/gm . 

DW' ® substrate 

= 9.0 ctn/hr 
CR| = 10 hr"'' 

Uq = 10 or 20 depending upon initial DSD 

Ug = 15 


* -jf * 
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iiPPMIDIX F 

JOLGOBITIM mn MPHTE-CAELO SIMI7L ATI0H 

Ibpulation balance equation (E, l) when laslved with the 
growth equations (E.4- to E.7)» is capable of giving rise tD tie 
growth behavior in n-alkane fermentation systems. Bbwever, 
analytical solution to these equations is very complex due to 
the following considerations: 

1. i, Ibpulation balance equation is coupled with the 

growth equations. 

ii. The growth equations for the dispersed phase are 
coupled in themselves and with those for the conti- 
nuous phase. 

2. i. Population balance equation is a nonlinear integra- 

partial differential equation with the population 
density, f, appearing under multiple integral signs 
in two different terns with variable limits, 
ii. The growth equations for the dispersed and the 
continuous phases are also non linear. 

The prediction of growth behavior of gas-oil fermentation 
systems is one order of magnitude more complex becsuse the 
population density of drops in this case is a trivariate function, 
the variables being v (the volume), S ( the substrate concentration), 
and n (the surface biomass density). 
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Monte— Carlo siEUlation procedure can, liowever, be 
conveniently used to analyse the problem under consideration, 
provided care is taken to see that the simulation procedure 
reproduces the physical process as closely as possible. It can 
be (tone by following exactly the course of events as is proposed 
by the model* The method is not only simple but also provides 
for incorporation of more complicated models for the particulate 
phenomena, A Monte-Carlo simulation procedure suitable to the 
problem at hand was developed by Shah (1974) and will be utilized 
here. 

The Methodology of Simulation ; 

Hydrocarbon fermentation systems fall into a general 
category popularly known as ’particulate systems* , The systans 
consist of particles undergoing random phenomena which affect 
the properties of these particles. The random phenomena could be 
either monoparticle, i.e. those involving only one ^article at a 
time, o r mul ti particle, Ibcanples of mo no par tide events are 
breakage, birth and death, entry aid exit from the system, etc. 

Of multijarticle events, only biparticle random phenomena like 
two particle coalescence, collision, etc. are imxiortant. Events 
involving more than tvso particles at a time are not very frequent 
and notof much interest in systems of particles like drops and 
bubbles. The randan phenomena are specified in terms of frequencies 
of their occurrence which may be functions of properties of the 
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particles taking pert in the event, and of environmental 
variables like speed of agitation, power input, dispersed phase 
fraction, temperature, shape and size of the vessel, and the 
agitator, etc. The random phenome na may be accompanied by growth 
of particles also. Both of these affect the state of the system. 
Growth causes a continuous change in state while change due fo 
random phenomena occurs abruptly at discrete time points at 
different intervals. Thus, the specification of the dynamic 
behavior of the system requires statement of (a) initial conditions, 
( b) frequency of random events, and ( c) grDwth behavior of the 
particles. 

State of the particulate system at any time is specified 
by the number of particles and their individual properties. 

Evolution takes place by a sequence of random events between which 
growth takes place. The time between the two successive random 
events is a random variable and is termed 'interval of quiescence' 
(IQ). IQ is the period during which no random event takes place 
and its probability distribution depends upon the frequency of 
events and tlB properties of existing particles. The simulation 
proceeds by generating a sequence of IQ and updating the proper- 
ties of tte particles to the end of each IQ when a random event 
takes place and causes sn abrupt change in tte properties of the 
particles involved. It is assimed that the random event occurs 
infinitely fast. A aifficient number of repetitions will ' 
allow the estimation of the required quantities. 
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file hydrocarbon fermentation model under consideration 
envisages only co al e scenes- red isperaion events accompaiied by 
growth of microbes and consumption of substrate, These random 
events are biparticle in nature. The frequency of this event is 
denoted by q and may depend upon the properties of the tv/o 
particles involved. Thus, 

q = q (x(t) , y( t)) (F.1) 

It is defined by, 

Pr l^a given particle of property x interacts with another 
particle of pirperty y in the time interval t and t+dt 
= q.(x, y) dt + o(dt) (P.2) 

Let the state of the syston at any time be described as, 

|xj^(t), i = 1,2,...,ll(t)| (P.3) 

where N( t) is the total number of particles at time t, and x^ 
is the property of the ith particle, H^(t) remains constant in 
the present investigation due to the specific nature of the drop- 
interactions proposed, q is also assumed to be constant and sane 
for all the particles. The algoritlm for this simulation has been 
developed in detail by Shah (1974), 

The probability distribution of IQ may be derived in the 
following manner; 

Let T be a random variable representing time between two 
successive biparticle events. If P( is 1die probability 



173 


that no biparticle event occurs between time t and t+'V', given 
the state of the system at time t, 

nr/A^) = p(T^r/^) 

Prom the definition, 


nr* arj ) = P(t| )• (no biparticle event occurs 




= ny 


between 'J' and + d3^ ) 

lT(t)-1 l[(t) 


- E E 4 

3=1 k=3+1 ^ 


x.(t),x,(t)> d'J' 


(P.5) 


The summation term accounts for all the possible pairs of 
particle population and will be valid only if F(t)3^2. 
IJividing by d and taking the limit d'^-^O, 


L. 

dy 


nr 


lT(t) 


‘A. 


1 ( t )-1 

) =- P(Tj ) q.x^(t),Xj^(t) ■ 

: h 3=1 k^+1 ^ 


(P.6) 


which leads to. 


t+T , . / > 

f- ... N(t)-1 F(t) 


nji ) = exp 


A^ 


- / Y. H 


3=1 k=3+1 


I dt’ 
(P.7) 


because by definition, 


P(0 j ) = 1. 


The cumulative probability distribution of T is related to P by. 
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KTI ) = Pr (T ) 


1 - Pr 








1 - p 


^l, 

L "'-t. 


Therefore 


"fc+37i 


'l(t)-1 N(t) 

ny ) = 1-exp I- j Y Y 

- I L_ L_ 

t 3=1 k=3+1 




}«'] 


(P.8) 


Similarly, the conditional probability for a particular pair 
(ijj) to be involved in the biparticle event is obtained as 
fo Hows : 

Define an indicator variable, u, such that 

u = if tl® particles of properties and x^ undergo 

the bi particle event 

Since the required probability is proportional to q(x^, x-), 

?/e have , 


Pr u = (i, 3 ) 


- 

h rp J 


N(t)-i Nrt) 


Z_ 

k=1 


r 

l=k+1 


(P.9) 


The coalesced drops again redisperse into two daughter 
droplets whose sizes could be |x| (x^+x^-x)| . The pro babil ity 
for formation of a dau^ter droplet of size x is given by. 



175 


M 


Pr 


X< X 


A^,T,u = (i,3)J 


X 


Xi + X. 


(P.10) 


Startixig fxom an initial state, realization of tte 
population process of a batch two— phase microbial propagator 
consists of sampling of the following representative sequence 

..... 


* * * 
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